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摘要 : 1/1 共振 是 限制 性 三 体 问 题 或 行星 三 体 问题 中 最 低 阶 的 平 运动 共振 ， 平 面 圆 型 限制 性 三 
体 问 题 中 的 5 个 平 动 点 对 应 1/1 共振 的 5 个 平衡 点 ， 束 缚 在 不 同 的 1/1 共振 态 时 ， 系 统 在 空间 中 
表现 出 不 同 的 运动 轨迹 ， 如 晴 时 型 轨道 、 马 蹄 型 轨道 等 ， 而 在 大 偏心 率 和 高 倾角 下 ， 系 统 存在 其 
由 类 型 的 1/1 共振 态 。 太 阳 系 内 存在 许多 自然 共 轨 系统 ， 如 木星 的 特洛伊 小 天 体 及 土 卫 的 共 轨 
运动 等 ， 近 年 来 ， 在 近 地 空 间 中 己 探 测 到 数 十 颗 与 地 球 共 轨 的 小 天 体 ， 在 轨道 的 长 期 进 动 、 大 行 
星 摄 动 等 因素 的 影响 下 ， 这 些小 天 体 可 能 在 不 同 的 1/1 共振 态 之 间 频 繁 地 迁移 。 哈 密 顿 系统 和 
Hill 三 体 问题 是 研究 1/1 共振 的 两 种 有 效 途径 ， 按 照 一 般 摄 动 理论 的 常规 做 法 ， 用 哈密 顿 系统 处 
理 1/1 共振 的 主要 内 容 有 : 摄 动 函数 展开 、 正 则 变换 或 平均 化 消除 短 周 期 项 、 理 想 共振 模型 及 长 
期 摄 动 理论 等 。 在 合适 的 摄 动 函数 下 ， 哈 密 顿 系统 可 以 研究 1⁄1 共振 的 全 局 动力 学 ， 另 一 方面 ， 
Hill 问题 主要 处 理 两 个 小 质量 天 体 的 相对 距离 较 近 时 的 相对 运动 。 限 制 性 问题 中 的 Hill 方程 是 人 
门 熟知 的 ， 而 具有 完全 相同 形式 的 一 般 Hill 三 体 问 题 可 以 处 理 两 个 小 天 体质 量 均 不 为 0 的 情形 。 
在 一 定 的 假设 下 ， 一 般 Hill 问题 的 运动 方程 可 以 写成 与 哈密 顿 系统 的 正则 方程 上 相似 的 正则 形 
式 ， 由 此 可 借助 平均 化 方法 建立 长 期 摄 动 理论 。 除 分 析 理 论 之 外 ， 亦 可 采用 数值 途径 去 研究 1/1 
共振 ， 如 周期 轨道 族 的 数值 延 拓 、 庞 加 莱 截 面 以 及 平衡 点 不 变 流 形 的 计算 等 。 
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三 体 问 题 研究 的 是 3 个 质量 分 别 为 mo, m, mz 的 质点 在 相互 引力 作用 下 的 运动 。 
当 mo > m, m 时 ， 主 天 体 mo 5 m R m 之 间 的 引力 远大 于 m 与 ma 之 间 的 引 
力 ， 此 时 ，mi 与 m2 在 空间 中 的 运动 可 视 作 受 到 摄 动 的 开 普 勒 轨道 ， 故 mi 与 m2 的 


收 稿 日 期 : 2021-04-12; ” 修 回 日 期 : 2021-09-26 


资助 项 目 : 国家 自然 科学 基金 (11773017, 11973010) 


通讯 作者 : 侯 锡 云 ，houxiyun@nju.edu.cn 


ChinaXiv 合 作 期 刊 


1 期 Ey, |: 限制 性 及 行星 三 体 问 题 中 的 共 轨 运动 研究 回顾 59 


运动 可 从 摄 动 理论 出 发 作 定性 或 定量 的 研究 。 对 三 体系 统 ， 可 定义 与 质量 有 关 的 参数 
u = (mi + ma)/(mo + mi +m), p « 1 Ë. m ss 的 情形 被 称 作 行星 三 体 问题 ， 而 
当 m/m 趋 于 0 或 无 穷 的 情形 被 称 作 限制 性 三 体 问 题 。 在 太阳 系 中 ， 太 阳 (mo)- 小 天 体 
(m1)- 行 星 (mo) 组 成 的 三 体系 统 的 运动 可 近似 为 限制 性 问题 ， 太 阳 (mo)- 行 星 1 (mi)- 行 星 
2 (m3) 组 成 一 个 行星 三 体系 统 。 容 易 想 象 ， 限 制 性 或 行星 三 体 问题 模型 在 系 外 行星 、 系 乡 
小 行星 、 三 星 (triple stars) 等 自然 天 体系 统 的 研究 中 同样 可 以 找到 适用 的 场景 。 共 轨 运 动 
是 三 体系 统 中 一 类 特殊 的 运动 ， 其 特征 是 mi 与 ma 的 轨道 半 长 径 相 当 (a, © a2)， 若 mi 与 
m» 的 轨道 偏心 率 和 轨道 倾角 都 很 小 ，mi 与 mo 便 几 乎 在 同一 轨道 上 运动 ， 这 解释 了 “ 共 
轨 运 动 ”的 含义 。 对 共 轨 运动 更 精确 的 说 法 是 1/1 平 运动 共振 ， 即 mi 与 ma 的 平 经 度 之 差 
(Ai 一 A2) 在 某 个 值 附近 摆动 ， 此 时 Ay — Ao 的 一 阶 导 数 (ni — na， 即 平 运动 角速度 之 差 ) BE 
近 于 0， 和 且 mi 5 ms 的 轨道 周期 相当 ， 由 开 普 勒 第 三 定律 ， 这 对 应 着 前 面 提 到 的 共 轨 运动 
的 特征 ;因为 n1 : nz 污 1:1， 故 称 系统 此 时 处 在 1/1 平 运动 共振 ， 且 称 Ai — A2 为 1/1 Fiz 
动 共振 的 共振 角 。 


= 为 了 结合 实际 的 自然 天 体系 统 ， 人 们 一 般 关 注 的 是 限制 性 或 行星 三 体系 统 在 某 参 考 系 
e 下 稳定 的 共 轨 运动 ， 这 对 应 着 1/1 平 运 动 共振 中 的 稳定 平衡 点 ， 当 共振 角 Xi 一 AS 绕 着 不 
N 同 的 平衡 点 摆动 时 ， 系 统 在 物理 空间 中 表现 出 不 同类 型 的 轨迹 ， 根 据 共振 角 的 振幅 大 小 ， 
Cr 不 同类 型 的 运动 对 应 着 1/1 共振 局 部 或 全 局 的 动力 学 性 质 。 在 平面 圆 型 限制 性 三 体 问题 
= Hy m 在 一 个 固定 的 圆 轨 道上 绕 转 mo， 此 时 在 一 个 与 ms 共同 绕 转 的 旋转 坐标 系 下 ( 称 


为 会 合 坐 标 系 ) 研究 m 的 运动 是 合适 的 ; FERIAS P. mi 的 运动 存在 5 个 平衡 点 ， 
即 人 们 熟知 的 3 个 欧 拉 共 线 平 动 点 (Li, Lo, Ls) 及 2 个 拉 格 朗 日 三 角 平 动 点 (La, Ls), 24 
p < ue = 0.038 52 IJ, L, 和 Ls 在 线性 意义 上 是 稳定 的 ， 系 统 存 在 着 绕 转 Za (或 Ls) Big 
时 型 轨迹 以 及 同时 绕 转 Ls, La, Ls 的 马蹄 型 轨迹 ，7na 是 系统 的 一 个 物理 碰撞 奇 点 ， 而 事实 
E, p< 0.047 7 时 系统 存在 着 一 族 逆 行 绕 转 ma 的 稳定 周期 轨道 ”， 称 为 f 族 周期 轨道 ， 这 
一 族 周期 轨道 上 偏心 率 较 大 的 轨道 在 平均 化 系统 中 表现 为 位 于 Xi 一 à = 0 的 1/1 共振 稳定 
平衡 点 ”， 而 绕 转 这 一 稳定 平衡 点 的 轨道 常 称 为 拟 卫星 轨道 。 图 四 画 出 了 典型 的 蜗 虹 型 轨 
迹 、 马 蹄 型 轨迹 以 及 拟 卫星 轨迹 。 当 m 为 近 圆 轨道 且 小 天 体 m 为 近 圆 、 近 平面 轨道 时 ， 


= 系统 仍 存在 上 面 所 述 不 同类 型 的 共 轨 运动 ， 但 轨道 的 稳定 性 及 稳定 区 域 在 轨道 偏心 率 的 摄 

< 动 下 可 能 发 生变 化 。 而 当 m 的 轨道 偏心 率 和 倾角 较 大 时 ， 系 统 可 能 存在 不 同 于 以 上 类 型 的 

ww 1/1 共振 态 ， 例 如 绕 转 Ly, Ls 以 及 ma 的 大 振幅 1/1 共振 轨道 ”， 以 及 在 不 同 的 1/1 共振 态 
之 间 迁 移 的 跃迁 轨道 ”。 

行星 三 体 问题 较 限制 性 问题 因 系统 自由 度 增 加 而 更 为 复杂 。 行 星 三 体 问题 中 常用 的 坐 


标 系 有 两 个 : 日 心 坐标 系 和 雅 可 比 坐 标 系 ， 而 在 共 轨 运动 的 研究 中 常用 的 是 日 心 坐标 系 ”。 
值得 注意 的 是 :由 日 心 华 标 及 共 罗 动量 导出 的 正则 开 普 勒 根 数 下 ， 人 们 熟知 的 拉 格 朗 日 平 动 
点 不 再 对 应 着 圆 轨道 ， 而 是 有 一 定 偏心 率 的 进 动 梢 圆 ” 。 而 在 一 些 数值 研究 中 ””， 使 用 日 
心 坐标 及 日 心 速度 或 由 此 导出 的 开 普 勒 根 数 在 物理 上 更 加 直观 。 在 一 个 “合适 ”的 转动 坐标 
APA (可 以 是 z 轴 与 ma 固 连 的 质心 坐标 系 ”， 也 可 以 是 m 与 ma 之 间 的 平均 轨道 ” ), 
星 三 体 问 题 中 存在 着 与 限制 性 问题 类 似 的 共 轨 运动 类 型 。 如 在 近 圆 近 平面 情形 下 ， 拉 格 朗 日 


`N 
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图 1 


点 附近 存在 映 昱 型 轨道 和 绕 转 Ls, La, Ls 的 马蹄 型 轨道 ， 以 及 偏心 率 较 大 时 的 拟 卫 星 轨 道 。 


在 mi = mz, al = a», el = eə 情形 下 ，Giuppone 等 人 


注 : 红色 五 角 星 为 中 心 天 体 720， 蓝 色 五 角 星 为 行星 me. 


平面 圆 型 限制 性 三 体 问 题 中 量 时 型 轨道 、 马 蹄 型 轨道 及 拟 卫 星 轨 道 在 会 合 坐 标 系 X Y 平面 上 的 投影 


下 通过 数值 模拟 发 现 ， 在 三 角 平 动 点 


处 存在 一 族 反 拉 格 朗 日 周期 轨道 ， 以 L, 点 为 例 ， 该 族 轨 道上 m, ma 的 近 星 点 的 平 经 度 相 
差 为 240*， 而 不 再 是 60*。 与 限制 性 问题 类 似 ， 行 星 三 体 问题 中 1⁄1 共振 轨道 在 不 同系 统 参 
数 下 的 稳定 性 及 稳定 区 域 是 研究 关注 的 重点 之 一 。 
太阳 系 内 已 发 现 许多 1/1 共振 系统 。 表 四 列 出 了 与 大 行星 相关 的 共 轨 小 天 体 数 目 ， 其 中 
晴 时 型 轨道 的 数据 来 源 于 IAU Minor Planet Center (MPC)B, 而 处 在 其 他 类 型 轨道 的 小 天 


体 尚 未 得 到 IAU MPC 的 确认 ， 部 分 观测 相关 的 参考 文献 可 见 了 


回 
F Greenstreet 等 人 ”的 文章 。 


1/1 共振 系统 也 存在 于 太阳 系 行星 的 卫星 系统 中 ， 土 星 的 卫星 中 目前 已 知 存在 五 个 1/135 
die AZ: Calypso 和 Telesto 相对 于 Tethys 的 运动 为 赂 时 型 轨道 ，Helene 和 Polydeuces fH 
对 于 Dione 的 运动 也 为 赂 时 型 轨道 ，Janus y 与 Epimetheus 的 相对 运动 为 马蹄 型 轨道 (交换 
轨道 )。Laughlin 和 Chambers 首先 研究 了 在 太阳 系 外 发 现 共 轨 行 星 的 可 能 性 ， 并 研究 了 
稳定 的 共 轨 构 型 对 应 的 系统 参数 范围 。Giuppone 等 人 "研究 了 通过 径 向 速度 法 探测 太阳 系 
外 共 轨 行星 的 可 能 性 。Gozdziewski 和 Konachi 通过 分 析 径 向 速度 信号 指出 HD128311 和 
HD82943 可 能 是 行星 共 轨 系统 的 两 个 宿主 ， 但 这 一 可 能 性 未 得 到 后 续 观 测 的 证 实 ”。 事 实 


上 ， 


ELE Hn, IAE IRE SEVA BEER TA (包括 开 普 勒 空间 望远镜 


= 还 未 在 太阳 系 外 确认 


°nttps : //minorplanetcenter.net/iau/lists/Trojans . html (accessed on 2021 June 14) 
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发 现 共 轨 行星 系统 T, 


表 1 目前 已 知 的 与 太阳 系 大 天 体 相 关 的 1/1 共振 系统 数目 


行星 晴 时 型 轨道 ”马蹄 型 轨道 ” 拟 卫星 轨道 跃迁 轨道 
金星 一 一 1 1 
地 球 1 5 5 4 
火星 9 = - _ 
木星 > 10 000 一 = 1 
ae 一 一 一 4 
RE E 1 = 一 = 
海王 星 28 一 = = 


TE: 跃迁 轨道 指 在 不 同 的 1/1 共振 态 间 之 间 迁 移 的 轨道 。 


三 体 问 题 可 以 从 哈密 顿 系统 的 角度 来 研究 ， 该 方法 能 够 比较 容易 地 研究 1/1 共振 全 局 或 
局 部 的 动力 学 性 质 。 主 要 过 程 是 通过 正则 变化 尽 可 能 地 降低 系统 的 自由 度 ， 同 时 使 得 变换 后 
的 系统 仍 保留 共 轨 运动 主要 的 动力 学 性 质 。 但 是 ， 哈 密 顿 系统 仍 有 一 些 技 术 性 难题 需要 解 
决 ， 如 靠近 碰撞 奇 点 时 展开 式 的 收敛 性 问题 、 共 轨 运 动 的 全 局 角 作用 量变 量 的 显 式 表达 式 
等 。 研 究 共 轨 运 动 男 一 有 效 的 方法 是 Hill 的 三 体 问 题 模 型 ， 在 这 一 模型 下 ， 共 轨 运 动 可 分 
为 两 部 分 考虑 : (1) 外 部 近似 ， 当 m, 与 ma 距离 较 远 时 ， 系 统 由 两 个 独立 的 开 普 勒 轨道 近 
似 ; (2) 内 部 近似 ，mi 与 ma 距离 较 近 时 ， 系 统 的 运动 由 Hill 方程 描述 。 除 了 以 上 两 种 模型 
为 代表 的 分 析 方 法 ， 数 值 积分 及 其 衍生 的 结果 是 研究 共 轨 运动 最 直接 的 途径 ， 也 是 检验 不 
同 理论 的 基石 ， 基 于 数值 积分 的 周期 轨道 延 拓 和 庞 加 莱 截 面 是 研究 相 空间 结构 的 有 力 工具 。 
第 回 章 对 哈密 顿 系 统 、Hill 三 体 问 题 等 共 轨 运动 研究 中 常用 方法 作 简 要 介绍 ， 第 回 章 对 共 轨 
运动 相关 的 动力 学 研究 作 简要 回顾 ; 第 田 章 总 结 本 文 介绍 的 内 容 ， 并 对 未 来 可 能 的 研究 方向 
作 简 单 的 讨论 。 


2.1 REMAR 
2.1.1 "à E WB S AK 

(1) 限制 性 三 体 问 题 

id ri, T2 分别 为 m1, ma 相对 于 mo 的 日 心 距离 矢量 ，~ 为 m 的 模 。 在 一 个 XY 平面 
与 ma 轨道 平面 重合 的 惯性 坐标 系 下 ， 由 牛顿 万 有 引力 定律 ， 不 难 推导 出 小 天 体 m 的 运动 
方程 如 下 : 


T2 T; — T° 
= T1 : Gm» & 1 
dt? ri | ( r3 71 一 =) ü) 
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由 式 (1) 容易 写 出 系统 的 拉 格 朗 日 函数 : 
ee ee (5 "E : ) . (9 
2 Ti T5 Iri — rə| 
H dv BJ H eG cn] de S: Hi # 2 BJ Hz] E A RED EL s 
OL . 
r = Əri —maT; , (3) 
Hi = Pr $ Ti L = Lou Gms (mo Eu ma) | Gm m2 72 71 1 s; (4) 
1 2 rı rj [rı 一 7a| 


考虑 到 m, = 0， 对 哈密 顿 系 统 作 以 下 等 价 变换 : 


Hı > Hi/mi , p,, > pr /mI 


则 得 到 限制 性 三 体 问题 的 哈密 顿 函 数 ; 


1 1 
Hy = Pa - E + Gm [Ir eos - X) ， (5) 
其 中 ， 
Ti ' To 
Pri =i, 1 =G(Mmo +m), cosy = a A=]|ri-r|. (6) 
1/2 


由 小 天 体 m1 BUSES As SITE ri, pu, PHS BITS ARE (a, e, i, 2, w，M)， 这 一 变量 
变换 不 是 正则 的 ， ri, p,, 与 Delaunay 根 数 之 间 的 变换 是 正则 的 ，Delaunay ff 


aa Gi =l /l—e2, Hi = Gi cosi (Y) 
h=M, gw, h = fA 
以 上 变量 在 偏心 率 为 零 时 不 是 良 定 的 ， 为 此 我 们 采用 无 奇 点 的 形式 : 
ame Ch = B, (yi-4-1) gH TEE, m 
MA = Mice fA, g =w tR, hi — fA 
由 于 行星 ma 在 固定 的 开 普 勒 椭 贺 上 绕 转 mo， 类 似 地 可 给 出 平面 无 奇 点 Delaunay 根 数 : 
La = V2a2 = const , G5 = La (v3 一 1) = const (9) 


入 2 一 入 2.0 十 Tot " gu = const 


HP, po = G(mo + ma), nə = (ua/a3)0/2. sÑ (5) 通过 ro 显 含 时 间 t， 可 人 为 地 引入 一 个 
新 的 自由 度 (Ao, 42) 得 到 一 个 自 洽 的 哈密 顿 系 统 : 


My 1 
Heer ud = Eu 1 
1 2a, + no À> + Gm» (5 cos Y% 5) (10) 
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(2) 行星 三 体 问题 

下 面 给 出 一 般 NN 体 问题 在 正则 日 心 4 
对 应 NN = 3 的 情形 。 在 一 个 空间 惯性 坐标 系 下 ， 记 (ui, ü, = m 
< 斩 动 量 ， 则 容易 写 出 系统 的 哈密 顿 函 数 : 


举 标 下 的 哈密 顿 函数 的 简单 推导 ”， 行 星 三 体 问 题 
Ui 


cu Am; 的 位 置 矢量 和 


= 


H,=T+V, T= pes s | (11) 


N-1 | (12) 
i=0 


系统 总 线 动 量 守 恒 ， To AEE, WADE T To = 0. 容易 验证 从 (ui, à) 到 (ri, Ti) 是 一 个 
FE 则 变换 ， 则 哈密 顿 函数 可 重 写 为 : 


E: I 


"n N-1 DH Gmom;j n Tic Pj i (13) 
< > 20; |r| > 7720 r;| š 
KH, ñ; = momi/(mo 十 mi;)。 容 易 看 出 ， 式 (03) 中 第 一 个 求 和 号 中 的 项 对 应 着 m; 的 开 普 


勒 运动 而 第 二 个 求 和 号 中 的 项 为 摄 动 项 。 记 (ai, €i; lis (2;, Wi, Mi) 是 1 (Ti, Til Bi) 决定 
的 开 普 勒 根 数 ， 则 哈密 顿 函数 可 重 写 为 : 


Hibi T;;T; Gmjm; 
H, = | š 14 
< = 2a; y» * Trio Ir; = Tj| ( ) 


值得 注意 的 是 ， 由 式 (12) 可 知 v; A m » MH (ri, 7i/ Bi) 决定 的 开 普 勒 椭圆 与 实际 速 
方向 并 不 相 切 ”。 与 式 (7) 类 似 ， 可 定义 行星 三 体 问题 页 下 的 Delaunay 根 数 : 


uae G,= L /1-e2, H,= Gicosi T 


l;—Mi, g=wi, hi = Ri 


为 解决 偏心 率 (轨道 倾角 ) 为 零 时 角度 量 wQ) 不 是 良 定 的 问题 ， 可 类 比 式 (8), (9) 对 幅 角 
进行 线性 组 合 ， 引 入 无 奇 点 形式 的 Delaunay 根 数 。 但 在 e = 0 或 i= 0 时 ， 无 奇 点 形式 的 
Delaunay 根 数 中 的 作用 量变 量 也 等 于 0， 运 动 方程 中 会 出 现 小 分 母 现象 ， 为 解决 这 一 问题 
可 引入 庞 加 莱 无 奇 点 变量 ， 其 具体 形式 不 再 叙述 
21.2 I BRR 

应 用 摄 动 理论 的 第 一 步 是 摄 动 函数 展开 ， 即 将 哈密 顿 函 数 中 的 摄 动 项 展开 为 开 普 勒 根 
数 或 其 他 正则 变量 的 容 级 数 形式 。 引 言 中 己 经 提 到 ， 共 轨 运 动 系 统一 般 满足 ar s a3) F 
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虑 到 较 小 轨道 偏心 率 ， 进 而 有 ri moros BEH my/ra (BR ro/r1) 作为 小 量 进行 展开 (B AG 
地 球 卫星 或 等 级 式 三 体 问题 中 的 Kaula 2URJEU 7) 一 般 并 不 适用 。 对 于 平 运动 共振 (包括 
1/1 共振 )， 经 典 的 展开 方法 是 考虑 系统 近 圆 近 平面 运动 ， 使 用 经 典 的 椭圆 展开 式 四 ， 首 先 
将 摄 动 项 展开 为 偏心 率 和 轨道 倾角 的 寡 级 数 呈 县。 由 限制 性 或 行星 三 体 问题 的 哈密 顿 函数 式 
(10) 和 (14) 可 知 ， 其 摄 动 项 可 分 为 两 部 分 : 直接 项 与 间接 项 ， 其 中 直接 项 即 为 包含 相对 距 
离 倒数 1/|r1 — ra| 的 项 。 事 实 上 ， 使 用 椭圆 展开 中 基本 关系 式 以 及 位 置 速度 矢量 与 轨道 根 
数 的 转换 关系 中， 对 间接 项 进行 展开 并 不 困难 ， 下 面 考虑 直接 项 的 展开 ， 类 似 于 式 (6)， 重 
新 记 直 接 项 为 : 


1 E 


a x i +r? 一 2rirjcos0) ° , (16) 
其 中 ， cos0 = TA, 假设 O(e;) = Ole) = Oli) = O(ij) 2 e < 1, MR (16) 可 表示 为 : 
il j 
=a pa 
Ai; B | e ^ 


ili 2 
A= a; + a; = 2a,a; COS Tjj 


4 =r; -a +ri— as (17) 
Zə = 2ajaj COS Tjj — 2riT; cos 0 
Tij = À; == À; 
由 假设 可 知 ，O(21) = 0(22) = e， 若 进一步 假设 O(A)= 1， 则 由 泰勒 展开 可 得 : 
一 工 = n 2n — 1)!! _ 2n+1 u n MS 
(A+ Z, + Z2) 2 = >` (—1) ( om A ° » BE Z . (18) 
n=0 1 一 0 


利用 几何 关系 式 及 椭圆 展开 式 ， 可 将 ZI ZU 展开 为 e; ey, d i; IRE KERR 
数 的 一 般 形 式 的 推导 ， 附 录 四 给 出 了 限制 性 三 体 问 题 及 行星 三 体 问 题 中 直接 项 与 间接 项 截 
WME Ole) 的 展开 式 。 在 此 种 展开 方式 中 ， 我 们 假设 a? + a2 — 2aiaa cos (Ay — A9) = O(1), 
而 这 一 条 件 对 于 共 轨 运动 并 不 总 是 满足 的 。 事 实 上 ，o2 + a3 — 2aja2 cos (Ay — A9) = 0 是 
(A.3) sti — 个 奇 点 ， 但 这 一 奇 点 仅 在 m 及 m 的 轨道 偏心 率 及 轨道 倾角 均 为 0 时 才 对 应 
着 系统 的 碰撞 奇 点 ， 因 此 ，az + a2 一 2aias cos (A1 — A2) = 0 是 经 典 展开 式 (A.3) 存在 的 一 
个 技术 性 问题 。 已 有 研究 指出 全 县， 受 该 奇 点 的 影响 ， 经 典 展 开 式 在 共 轨 运动 中 的 拟 卫 星 
区 域 失 效 。 式 (A.3) 另 一 个 固有 的 限制 是 偏心 率 须 小 于 拉 普 拉 斯 极限 ， 且 当 轨 道 倾角 较 大 时 
展开 式 收敛 效果 较 差 。 

R (A.3) 本 质 上 是 以 零 偏 心率 及 零 轨道 倾角 (e = i = 0) 为 参考 值 对 1/A 进行 展开 。 
Ferraz-Mello 等 人 在 对 太阳 系 小 天 体 平 运动 共振 的 研究 中 使 用 了 另 一 种 展开 方式 ， 即 使 用 
非 零 偏心 率 及 轨道 倾角 (e = eo, i = do) 为 参考 值 对 1/A 进行 展开 ， 在 平面 限制 性 三 体 问题 
F; Ferraz-Mello"" 针对 一 阶 平 运 动 共振 给 出 了 这 一 展开 方法 的 具体 算法 。Ferraz-Mello 和 
Sato" 将 展开 方法 推广 至 更 高 阶 平 运动 共振 ，Roig 等 人 四 进一步 推广 至 空间 情形 。 这 一 方 
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法 的 限制 在 于 : 展开 方法 本 质 上 是 局 部 的 ， 故 可 能 不 适用 于 较 大 振幅 的 共振 轨道 ”， 此 外 ， 
其 平均 化 需 借 助 数 值 平均 ( 见 下 一 节 )， 无 法 分 析 地 完成 。 
关于 摄 动 函数 直接 项 展开 还 需 提 到 的 一 点 是 ， 式 (A.3) 并 不 是 一 般 意义 上 的 完整 摄 动 函 

数 展开 ， 因 分 母 中 还 含有 cos (和 1 一 A2) 相关 的 项 。 然 而 ， 如 果 我 们 的 研究 目标 是 1/1 平 运动 
Life (HII Ay — Ap 是 共振 角 )， 式 (A.3) 已 经 足够 用 来 进行 平均 化 以 消除 快 变 量 ” ( 见 下 一 节 ); 
另 一 方面 ， 如 果 我 们 的 研究 目标 是 一 般 的 平 运动 共振 (BD pA1 — ghz 十 .… 是 共振 角 )， 则 还 
需 进 一 步 引 入 拉 普 拉 斯 系数 将 分 母 n 展开 为 傅 里 叶 级 数 形式 ”"。 还 需 提 到 的 一 点 是 ， 如 果 
我 们 关注 的 是 局 部 的 1/1 共振 态 ， 那 么 可 以 对 分 母 L, 在 共振 角 Xi — Ap 的 某 个 共振 中 心 处 展 
开 ， 此 时 容易 得 到 1/1 共振 局 部 的 解 ™。 
21.3 正规 化 与 平均 化 

引言 中 已 经 提 到 ， 处 理 哈 密 顿 系统 的 基本 思想 是 通过 一 系列 变换 (正则 变换 等 ) 尽 可 能 
氏 系统 的 自由 度 ， 同 时 使 得 变换 后 的 系统 保留 原 系统 中 我 们 关注 的 动力 学 性 质 。 严 格 地 
， 想 要 降低 系统 的 自由 度 ， 则 必须 找到 相应 的 运动 积分 ， 对 于 限制 性 及 行星 三 体 问题 则 比 
较 困难 。 引 言 中 已 经 提 到 ， 这 里 所 关注 的 共 轨 运动 本 质 上 是 一 个 近 可 积 系 统 ， 即 m, 与 ms 
绕 转 mo 的 开 普 勒 运动 受到 相互 引力 的 摄 动 。 摄 动 理论 中 处 理 这 样 的 系统 有 严格 的 方法 ， 即 
借助 哈密 顿 函 数 的 震级 数 形 式 ， 通 过 一 系列 正则 变换 将 原 哈 密 顿 函数 变换 为 一 个 截断 至 革 
个 阶 数 的 正规 化 系统 ， 同 时 保留 一 个 摄 动量 级 为 N + 1 的 余 项 。 若 以 上 变换 在 N > oo 
时 仍 收敛 ， 则 系统 自由 度 的 减少 是 严格 的 ， 若 以 上 变换 截断 至 某 个 有 限 阶 数 N， 则 变换 后 
的 正规 化 系统 可 视 作 原 系统 的 一 个 精确 至 O(eN) 的 近似 。 

与 正规 化 对 应 的 ， 降 低 系统 自由 度 的 另 一 方法 是 对 哈密 顿 函数 或 运动 方程 直接 进行 “ 平 
均 化 ”。“ 平 均 化 ”要 求 系统 不 同 的 自由 度 存 在 等 级 式 的 结构 ， 即 不 同 自 由 度 的 频率 (或 特征 
时 标 ) 的 大 小 在 不 同 的 量 级 ， 容 易 理 解 ， 在 较 长 时 标 下 看 ， 周 期 较 短 的 自由 度 可 视 作 短 周 其 
的 “噪声 ”， 对 系统 长 期 演化 贡献 的 “均值 ”为 零 ， 这 解释 了 “平均 化 ”的 含义 。“ 平 均 化 ” 
的 具体 做 法 是 对 哈密 顿 函数 或 运动 方程 中 短 周期 的 幅 角 进行 积分 平均 ， 变 换 后 的 系统 中 不 
再 出 现 短 周期 幅 角 ， 故 系统 自由 度 自 然 减少 。 需 要 指出 的 是 ， 对 于 多 自由 度 系统 ， 与 正规 化 


X dg 


( 方式 不 同 ， 平 均 化 系统 与 原 系统 之 间 的 关系 尚 无 严格 的 理论 ， 从 这 个 角度 看 ， 平 均 化 是 一 种 
= 经 验 性 的 做 法 ” 。 
PT 束缚 于 平 运动 共振 的 系统 有 一 个 共同 的 特征 :系统 不 同 的 自由 度 存在 等 级 式 的 结构 ， 即 


不 同 自由 度 的 频率 (或 特征 时 标 ) 的 大 小 在 不 同 的 量 级 。 在 1/1 共振 中 ， 共 轨 天 体 的 平 运动 
是 系统 中 周期 最 短 的 自由 度 ， 哈 密 顿 函数 中 的 开 普 勒 项 是 平 运动 的 可 积 近似 。 哈 密 顿 系统 的 
第 一 步 便 是 通过 正则 变换 或 平均 化 消除 平 运动 对 应 的 自由 度 。 将 附录 四 中 的 展开 式 代入 式 
(10) 或 式 (14)， 即 得 到 限制 性 或 行星 三 体 问题 的 哈密 顿 函 数 的 展开 形式 (不 妨 统一 记 为 H), 
为 引入 1/1 共振 角 ， 作 以 下 正则 变换 : 


h: (Li, `, Ao, A2) > (Li, T = A1 — Ag, J2 = (Li + 42), A2) ， (19) 


HH, Ap 对 应 着 系统 平 运动 的 自由 度 。 为 消除 这 一 自由 度 ， 可 应 用 Lie 变换 ”对 系统 进行 


N 
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正规 化 ， 记 正规 化 系统 为 8*， 由 生成 函数 可 得 到 映射 : 


F:a*>0o, (20) 


其 中 ，o 代表 旧 变 量 ，o* 代表 新 变量 。 由 式 (20) 可 得 : 


H*=Ho F, (21) 


H H* 中 不 再 显 含 Aš. AH Mathematica 等 符号 计算 软件 以 及 基于 Lie 变换 的 Deprit 3518 
算法 一 ， 在 计算 机 上 实现 哈密 顿 系统 正规 化 的 计算 并 不 困难 。 关 于 正规 化 最 后 需要 指出 
的 一 点 是 ， 在 对 哈密 顿 系统 H 正规 化 之 前 ， 需 首先 对 H 中 的 项 按 其 量 级 进行 整理 ， 若 考虑 
小 量 为 O(m,) = O(e;) = O(e;) = Ol) = O(ij) =e 之 1， 则 行星 三 体 问 题 的 哈密 顿 函 数 按 
量 级 表示 为 : 


H = Hj, + H> + Hgts: € (22) 


在 研究 中 ， 男 一 常用 的 做 法 是 ， 在 整理 万 中 的 项 时 仅 考 虑 质量 为 小 量 ， 即 O(m;) = e < 1, 
则 行星 三 体 问 题 的 哈密 顿 函数 表示 为 : 


H; = H3, T H3; š (23) 


容易 理解 ， 第 二 种 方法 对 量 级 的 整理 并 不 严格 ， 但 在 偏心 率 及 轨道 倾角 较 小 时 ， 式 (23) 会 
给 计算 带 来 方便 。 事实 上 ， 无 论 使 用 式 (22) 或 式 (23)， 其 正规 化 结果 均 含 有 无 穷 项 ， 记 式 
(23) 的 正规 化 结果 为 : 


H3” HS + Hy + HS 4€ . (24) 


juni 


Hak (20)， 记 由 生成 函数 得 到 的 新 旧 变 量 之 间 的 映射 为 F?。 
前 面 提 到 ， 对 哈密 顿 函数 H 直接 进行 积分 平均 可 以 得 到 平均 化 哈密 顿 H, 


= 1 _ 
H= = f Hoh 'dàs . (25) 


不 难看 出 ， 式 (25) 等 价 于 引入 共振 角 后 扔 掉 摄 动 函 数 表 达 式 中 显 含 Aə 的 项 ， 故 从 附录 四 推 
得 限制 性 或 行星 三 体 问 题 的 平均 化 哈密 顿 是 十 分 直接 的 ， 以 平面 行星 三 体 问 题 为 例 ， 其 平均 
化 哈密 顿 函数 为 ; 


H, = — me | Gmim2 (TN = T>) ; (26) 
AN rh, 
mMo 1 et 3E e ü 
Y= 1 27 一 人 2 
V (mo + m.) (mo + m3) J/a182 ( 2 EUNT aaa eel Pp PR 
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3a? 1 3a4 15 
T5 =l l e? 区 | ui 5203008 7 | a( i | a2a; 


3 9 
1 1 8 a50jcos T — 30207008 2r) + 


2 


3a2 1 Jas, 15 3 9 
e2 k: (32-5 cos) we ( a | 3 tae 5143 cos T— s dcos 25 


1 9 
" EZ COS (27— Ag) +7 ara2 cos agl + E (28) 


9 .. 8 _ 
€1€2 E (axa 4-aja) cos At (a2aj 4-aa1) cos (2r — Aj) 


h 


3 27 3 
+3430, cos (3r — Aj) J-a2a2 ( g Cos (r+ Ag) 本 cos (r Aa) 


其 中 ，A5 = 六 一 部。 前 面 提 到 ， 相 比 于 正规 化 ， 平 均 化 能 够 十 分 简单 地 降低 系统 自由 
无 法 严格 给 出 平均 化 系统 与 原 系统 的 关系 ; 尽管 如 此 ， 根 据 式 (24)， 不 难 发 现 : 


d 
m 


H; = Hy T H> . (29) 


即 在 最 低 阶 的 情形 下 ， 平 均 化 过 程 给 出 的 平均 哈密 顿 函数 与 正规 化 过 程 给 出 的 哈密 顿 函 数 
相同 。 式 (29) 一 定 程度 上 说 明了 平均 化 的 有 效 性 ; 除 此 之 外 ， 式 (29) 表明 在 一 定 精度 下 ， 
映射 r^ 建立 了 H, 与 平均 化 系统 H, 之 间 的 联系 。 图 回 给 出 了 平面 行星 三 体 问题 系统 平均 
化 的 一 个 算 例 : 给 定 初 始 平均 根 数 os. HX (20) 得 到 初始 瞬时 根 数 ogo， 对 完整 系统 H, 以 
Oo 为 初始 条 件 进 行 数值 积分 得 到 o(t)， 其 结果 对 应 图 Ba) 中 的 运动 轨迹 及 b) 一 d) 中 的 灰色 
粗 实 线 ， 对 平均 化 系统 He 以 os 为 初始 条 件 进行 数值 积分 得 到 二 (t)， 其 结果 对 应 b) —d) 中 
的 红色 虚线 ， 由 P) 及 平均 根 数 5(t) 可 得 到 F ort) HARIS bd) 中 的 蓝 色 或 黄色 
曲线 。 

关于 平均 化 还 需 提 及 的 一 点 是 ， 式 (25) 的 积分 平均 可 以 不 借助 摄 动 函 数 的 容 级 数 形式 ， 
而 是 直接 使 用 摄 动 函 数 原来 的 定义 ， 此 时 积分 平均 转 为 歼 曼 求 和 。 以 圆 型 限制 性 三 体 问题 为 
例 ， 式 (25) 转换 为 : 


N 
I 1 i ub = n 
M = jor s Pm yO Hoh 'o (e= 27) | (30) 


容易 理解 ， 在 式 (30) P, BR f A. 在 计算 中 从 0 到 2m 遍历 ， 其 他 轨道 根 数 ({o1,o2}/ 和 2) 
均 取 为 定 值 ， 故 当 N 取 充 分 大 的 有 限 值 时 ， 式 (30) 右 端 给 出 了 平均 化 哈密 顿 函 数 
Hio, 62}/d2=const 的 近似 值 。 通 过 遍历 ai 及 7， 便 可 在 1/1 共振 项 图 上 画 出 平均 化 哈密 
函数 Hio, vs=const 的 等 高 图 实际 上 ， 这 样 得 到 的 等 高 图 能 够 反映 系统 相 空 间 不 
区 域 处 1/1 共振 的 性 态 。 基 于 以 上 方法 ， 图 回 给 出 了 圆 型 限制 性 三 体 问 题 下 的 计算 结 
» E] Ba) 给 出 了 平面 情形 下 el = 0.25 时 的 1/1 共振 相 图 ， 图 Bb) 给 出 了 平面 情形 下 
el = 0.3, i1 = 36.8° 时 的 1/1 共振 相 图 。 需 要 说 明 的 是 ， 用 分 析 平 均 法 得 到 的 平均 哈密 顿 函 
AX (例如 式 (26)) 去 画 相 同 的 1/1 共振 相 图 ， 计 算 速 度 相 比 于 数值 平均 方法 会 显著 提高 ， 然 
而 ， 在 靠近 伪 奇 点 at + a2 一 2aiaa cos (和 1 — Az) = 0 或 磁 撞 奇 点 的 区 域 ， 平 均 哈 密 顿 函数 给 


xm oup s= 
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200 
t 
c) 


HS 系统 参数 为 mı = 1074, ma = 1.2 x 10 1, a) 在 一 个 随 mi FERN BU bka) 2 £: 4 


HERR Fo mo 处 在 (0,0) 


CHEFE) my 处 在 (1,0) (红色 十 字 )， 黑 is 出 线 给 出 了 ma 相对 于 mi 在 XY 平面 内 的 运动 轨迹 ，b) 一 d): 


ma 的 正则 开 普 勒 根 数 随时 间 的 演化 ， 其 中 灰色 曲线 为 完整 系统 数值 积分 o (t), EERI 


e(t), b) 中 蓝 色 曲线 为 F° o F(t), c) 与 d) pes (黄色 ) 曲线 分 别 为 ma (ma) 
2 平面 行星 三 体 问题 系统 平均 化 算 例 


e=0.25 i-0* 


q 1.00 


e =0.3 


F 均 化 系统 数值 积 
的 轨道 根 数 F^ o F(t). 


i=36.8° 


0.968 180 3060 0960 
7/(°) 
a) 
注 ，wi = Q, = 0. a) 深 色 部 分 对 应 着 磁 撞 曲线 ， 即 my 与 ms 的 物理 碰撞 奇 点 ，r = 0 处 平衡 点 对 应 着 完整 
系统 下 的 六 族 周期 轨道 ， 平 衡 点 附近 运动 即 为 拟 卫星 轨道 忆 四 ，b) 在 较 大 轨道 倾角 处 i = 36.8"， 系 统 出 现 了 
绕 转 La, Ls 及 mi 的 1/1 共振 轨道 回 
3 ” 圆 型 限制 性 三 体 问题 数值 平均 得 到 的 1/1 共振 相 图 
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出 的 相 图 结构 与 图 四 会 有 定性 的 差异 ， 这 表明 现 有 的 分 析 方 法 在 奇 点 附近 的 适用 性 受到 限 
制 ， 仍 需 改 善 。 
2.1.4 理想 共振 模型 

在 消除 平 运动 自由 度 后 ，1/1 共振 角 对 应 的 自由 度 在 系统 中 周期 最 短 ， 若 要 通过 正则 变 
换 进一步 消除 这 一 自由 度 ， 首 先 需 在 哈密 顿 函数 中 提取 出 1/1 共振 的 可 积 近 似 ( 即 一 自由 度 
的 理想 共振 模型 ， 在 正则 摄 动 理论 ”中 常 被 称 作 Hori kernel)。 将 el = es = 0 RAR (26) — 
(28)， 行 星 三 体 问题 1/1 共振 的 一 个 可 积 近 似 (理想 共振 模型 ) 如 下 


MB Moo | G Mo cos(A; — Az) 
T TI TIl5 
204 2a2 V (mo +m) (Mmo +m) Vada 


Hİ = 


(31) 


1 
va? + a2 = 2a1a2 cos( 和 1 = =) 


心 


T Kı = iy Hı, Ko = D> / pa 
N 2J; = Ly — Lo = kı ya — ko yüz , 0 =A . (32) 
2J5 = Li + Ly = kı ya + Kova , 05 = M + À2 
容易 验证 (Li, Lo, À, A2) > (Ju, J2, 01,02) 是 一 个 正则 变换 ， 在 该 变换 下 ， 矶 仅 显 含 共振 
角 0; = Ai — X， 故 为 一 自由 度 可 积 哈 密 顿 系统 ，.]2 是 一 个 运动 积分 。 按 照 理 想 共振 模型 的 


处 理 方式 ”， 下 面 将 系统 在 共振 中 心 处 展开 ， 根 据 mi 及 ms 的 无 摄 开 普 勒 运动 ，01 的 共振 
中 心 定 义 为 : 
e T = Nə < 2 = à . (33) 
[S 下 面 计算 当 式 (33) 成 立时 J, 的 值 Je HR (32) 得 ; 
= 2 2 
fm a= (22) , a =(=) . (34) 
Ky K2 
`Ë 将 式 (34) 代入 式 (33)， 可 得 ; 
Jie =v, pa s . (35) 


定义 两 个 新 变量 J 与 如 下 : 


Jia =uJo+ J, u = (36) 


K1 + K2 


将 式 (36) 代入 式 (34) 得 : 


J J 
va = VG T os Li E — , (37) 
1 2 
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sei. mu Oy. (38) 
K1 K2 


由 式 (37) 不 难 验证 J = 0 正 是 由 式 (33) 定义 的 共振 中 心 


1 2 1 
a (Eze) - [op Ky Ka - (E) 
G2 工 一 Z ni my Ky ua) ` 
将 式 (37) 代入 式 (31), BIRTH atn cH 在 共振 中 心 7 — 0 处 展开 。 注 意 到 ， 此 时 
HÀ 右 端 有 3 个 小 参数 : 


M4 mM 

u, T aD 2 
由 计算 容易 发 现 ， 式 (31) 中 的 开 普 勒 项 代入 式 (37) 展开 后 ， 结 果 中 关于 wv 的 一 次 方 项 为 0; 
为 了 使 展开 后 系统 截断 至 最 低 阶 时 保留 原 系统 主要 动力 学 性 质 ， 令 Ou?) = O(o1) = O(o2) 


(这 一 技巧 源 于 Karl Bohlin” )。 我 们 略 去 展开 的 具体 细节 ， 给 出 H] 展开 至 最 低 阶 表达 式 ; 


2 
Hii = ; (G + zz) (kı | Ko) u? + hon 区 5 E =] : (39) 
我 们 目的 是 得 到 1/1 共振 可 积 近 似 的 角 作 用 量变 量 ， 进 而 通过 正则 变换 消去 这 一 自由 度 ， 
这 一 问题 现在 归结 为 求解 HIS 的 角 作 用 量变 量 。 然 而 ， 通 过 计算 不 难 发 现 ，HJi 的 哈密 
顿 - 雅 可 比方 程 的 解 无 法 简单 地 用 初等 函数 写 出 (有 观点 认为 通过 椭圆 函数 可 能 给 出 显 式 表 
达 式 “")， 这 是 目前 从 摄 动 理 论 出 发 分 析 地 研究 1/1 共振 全 局 动力 学 的 一 个 问题 ， 事 实 
上 ， 即 使 能 够 写 出 可 积 近 似 的 角 作 用 量 ， 因 摄 动 函数 中 1/1 共振 角 以 ,的 形式 出 现 ， 而 
不 是 一 般 的 傅 里 叶 级 数 形式 ， 使 用 正则 变换 将 系统 正规 化 到 高 阶 将 十 分 困难 。 一 个 例子 是 
Garfinkel” 的 广义 理想 共振 模型 理论 ， 在 该 理论 中 Garfinkel 使 用 超 椭圆 积分 对 一 般 的 一 自 
由 度 共振 问题 给 出 了 二 阶 解 ， 基 于 该 解 ，Garfinkel 给 出 了 平面 圆 型 限制 性 问题 下 共 轨 运动 
极 坐 标的 二 阶 和 解 ””， 该 解 具有 全 局 的 适用 性 。 
最 后 ， 我 们 回 到 1/1 共振 可 积 近 似 HJ. Egit T HI] 的 典型 相 图 ， 以 及 与 相 流 对 应 
的 会 合 坐 标 系 下 不 同类 型 的 运动 轨迹 。 由 图 4 可知 ，Hi 保留 1/1 共振 近 圆 近 平面 时 主要 的 
动力 学 性 质 ， 即 与 5 个 平 动 点 对 应 的 5 个 平衡 点 ， 以 及 不 同类 型 的 1/1 共振 态 (RRL AY IL 
E 迹 、 马 蹄 型 轨迹 )。 

= 21.5 ”长 期 摄 动 理论 
车 通过 正则 变换 能 够 消除 原 系统 中 1/1 共振 角 对 应 的 自由 度 ， 则 变换 后 的 系统 给 出 了 
(Gi, Hi, ji, hi) 的 演化 ， 即 描述 了 共 轨 运动 的 长 期 演化 。 然 而 ，2.1.4 节 已 经 提 到 ， 由 于 1/1 
EJ B REGEL (或 称 作 理想 共振 模型 ) 的 复杂 性 ， 经 典 的 摄 动 理论 难以 消除 1/1 共振 的 自 
度 。 考 虑 到 关于 哈密 顿 系统 介绍 的 完整 性 ， 下 面 我 们 首先 对 限制 性 三 体 问题 的 3/1 平 运动 
k 振 的 长 期 演化 理论 -一 主要 是 Wisdom " X Henrard 的 工作 ， 作 简要 介绍 ， 进 而 对 1/1 
tk 振 长 期 演化 理论 作 简 单 讨 论 。 
Wisdom 等 人 使 用 经 典 的 摄 动 方法 处 理 了 平面 椭圆 型 限制 性 三 体 问题 3/1 平 运动 共振 
的 长 期 演化 。 下 面 简单 介绍 其 基本 思想 和 方法 ， 记 3/1 平 运动 共振 的 自由 度 为 (5, ó = 


DS 


[t 
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1.5 
0.02 
1.0; 
0.01; 
0.5: 
B 0 > 0 
-0.01 -0.5r 
-0.02 -1.0r 
-150 -50 — 850 150 5 oA oC 0 | 1 
A,-A,/(°) x 
a) b) 


3E: 系统 参数 为 mo = 1, mı = 0.004, ma = 0.001. a) 绿色 曲线 是 从 点 Ls 生发 的 分 界线 (seperatrix), HU 
以 点 Ls 为 极限 的 同 宿 轨 道 ; 红色 ( 蓝 色 ) 曲线 为 从 点 Li (点 La) 生发 的 分 界线 。b) 中 心 天 体 mo 位 于 (0, 0) 
Ab, m, 位 于 (1, 0) 4b, XY 平面 上 的 轨迹 给 出 了 me 在 mo 和 m, 的 会 合 坐 标 系 中 的 运动 。 绿 色 曲 线 包 络 的 
区 域 为 绕 转 点 La 和 Ls 的 晴 旦 型 轨迹 ， 在 绿色 曲线 和 红色 曲线 中 间 的 轨迹 对 应 马蹄 型 轨道 。 由 于 偏心 率 较 小 时 
所 卫星 轨道 并 不 出 现在 平均 化 系统 中 四 ， 图 中 没有 拟 卫 星 区 域 。 


图 4 行星 三 体 问题 1/1 理想 共振 模型 Ho 的 相 图 以 及 在 轨 线 在 会 合 坐标 系 的 投影 


À - 32), AS (ji, e) 有 关 的 自由 度 为 (z, 由， 系统 的 哈密 顿 函数 写 为 : 


2 
Hs = E — 36 + A (z, y, €2) + B (z, y; €2) cos (ó — P (z, y)) + C (z, y, e2) cos ó * (40) 


RF, A, B,C, P 均 为 已 知 函数 ， 有 具体 形式 可 参考 文献 Bm Wap, it O(z) = O(y) = 
O(e3) = e, O(9) =1， 则 有 O(4) = O(C) = 8, O(B) = < Wisdom X£—2V/4B X ||| , loll < 
$， 这 实际 上 是 平 运动 共振 常见 的 特征 ， 即 不 同 自由 度 的 频率 (周期 ) 大 小 在 不 同 的 量 级 。 
z, y 在 (D, o) 的 特征 时 标 下 变化 十 分 缓慢 ， 则 可 以 近似 将 变量 x, y “冻结 ”， 即 取 为 常数 ， 
此 时 系统 退化 为 一 自由 度 。 若 将 H, 中 的 无 摄 项 (BI —u?/28? — 39) 在 3/1 共振 平衡 点 处 展 
JF (代入 变换 下 = p. +u FREF, MELA TEAR) 并 省 略 高 阶 项 ， 则 系统 退化 为 一 自由 度 
E: 


Hj! = —ou? + B (z, y, eo) cos (ó — P (z, y)) , (41) 

其 中 ，a AUER. AAMA IIl 的 哈密 顿 - 雅 可 比方 程 可 写 出 通 解 ， 具 体 细节 可 见 

文献 BT, HH 引 。 简 言 之 ， 我 们 能 够 由 哈密 顿 - 雅 可 比方 程 得 到 生成 函数 S (0, 7)， 进 而 得 到 Hi 
的 角 作 用 量 (I, 9)， 以 及 3/1 共振 的 显 式 解 : 

P(t) = B(1, 0(t)), Pt) = elU, 0 的 ) . (42) 


BUA HÌ 给 出 的 3/1 共振 运动 与 完整 系统 Hs 的 结果 相差 不 大 ， 便 可 利用 式 (42) 对 哈密 
顿 函 数 Hs 或 由 Ha 导出 的 哈密 顿 正则 方程 进行 积分 平均 。 由 于 9 = X 一 3A 既 可 能 处 在 循 
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环 也 可 能 振动 ， 积 分 平均 此 时 写 为 : 
T Pmax 
_ 1 = de 
og [fto f f (43) 
0 pmin 
平均 化 哈密 顿 冰 数 为 : 
(H3) = E (z, Y, e2) (cos $) + F(a, Un e2) (sin 9) , (44) 
KB, RACE, F BEA, B, C, P EH. Wisdom 附录 中 用 正则 绝热 摄 动 理论 ”给 出 了 
平均 化 系统 (Ha) 更 严格 的 推导 ， 其 基本 思想 是 将 变量 z, y 的 恒 等 变 换 加 入 生成 函数 S, HU 
Si=S+ay. (45) 


L 


在 由 生成 函数 S 定义 的 正则 变换 下 ， 可 严格 地 导出 新 的 哈密 顿 函 数 HS, Wisdom 证 明了 当 


H3 截断 至 低 阶 时 (不 妨 记 为 HY 


整 系统 H; 的 一 个 “ 拟 积 分 ”)， 


)， 由 5S 定义 的 作用 量 工 恰好 是 系统 HS 的 运动 积分 ( 故 是 完 
HH H; 导出 的 哈密 顿 方程 与 平均 化 运动 方程 形式 上 完全 一 


致 ， 这 说 明了 在 哈密 顿 函数 或 哈密 顿 正则 方程 层面 进行 平均 化 是 合理 的 。 
前 文 已 经 提 和 到， 降低 系统 自由 度 必须 找到 相应 的 运动 积分 ， 对 于 近 可 积 的 系统 ， 则 需要 
找到 其 拟 积分 ， 在 上 面 介绍 的 Wisdom 等 人 的 工作 中 ， 作 用 量 了 是 原 系 统 3/1 共振 自由 度 对 


应 的 一 个 拟 积分 。 针 对 二 自由 度 的 近 可 积 系统 ，Henrard 等 人 在 经 典 的 摄 动 方法 基础 上 提出 


了 一 种 “半数 值 ”的 摄 动 理论 ” ,“ 半 数值 ”的 含义 是 指 计算 拟 积分 的 理论 基础 就 是 经 典 的 
摄 动 理 论 ， 但 Henrard 等 人 不 再 局 限于 得 到 拟 积分 的 分 析 表 达 式 ， 而 是 试图 对 原 系统 任意 
条 轨道 用 数值 方法 在 一 定 精 度 下 得 到 其 拟 积分 的 具体 数值 。 事 实 上 ，Henrard 等 人 给 出 的 


算法 能 够 计算 系统 


个 三 维 等 能 流 形 的 


个 横 截 面 上 任 一 点 的 拟 积分 数值 。 我 们 对 这 一 思 


路 作 简单 的 解释 : 二 自由 度 哈密 顿 系 统 若 存在 另 一 对 合 积分 (不 妨 也 记 为 门 ， 则 系统 可 积 ; 


一 般 情况 下 ， 可 积 系 统 的 相 流 微分 同 胚 于 二 维 环 下 


“， 考 虑 一 个 与 相 流 相 横 截 的 二 维 平面 ， 


则 微分 同 胚 于 三 维 环 面 的 相 流 在 二 维 截 面 上 的 投影 是 一 个 一 维 流 形 ， 而 不 同 轨道 (投影 的 一 
维 流 形 ) 对 应 的 积分 的 取 值 必 不 同 ， 故 积分 在 二 维 截面 上 的 等 高 线 即 给 出 了 等 能 面 上 的 


相 流 结构 ， 若 选取 合适 的 截面 ， 
的 长 期 演化 相 图 。 


H, = Ho (q 


其 中 ， (q, q2, Pi; p2) 是 正则 变 
初始 条 件 (quo; 92,0) P1,0; Da2.0)， 


的 周期 ， 则 系统 角 作用 量 分 别 为 


则 此 时 只 需 计 算 截 面 上 任意 一 点 积分 了 的 取 值 便 可 得 到 系统 


Henrard 等 人 考虑 的 二 自由 度 哈 密 顿 函数 形式 如 下 : 


1, Pi, M (q2, p2)) +EF (m, Di, q2, P2) - (46) 


量 ，s 是 小 参数 ， Ho(qi, Pi; M (q2, p»)) 可 分 离 故 可 积 ， 给 定 
系统 Ho 决定 了 一 条 ( 拟 ) 周期 轨道 ， 记 工 是 自由 度 (qi, pi) 
E 


t 

6, =2 
1 T, 
" 

0, = 2 
2 T, 
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般 而 言 ， 只 有 Ho 的 形式 非常 简单 时 (例如 谐振 子 )， 以 上 积分 才 可 能 显 式 地 解 出 。Henrard 
等 人 指出 ， 当 Ho 形式 十 分 复杂 时 ， 可 以 从 计算 的 角度 出 发 ， 利 用 式 (47) 计算 一 条 给 定 的 
(B) 周期 轨道 的 角 作用 量 的 具体 数值 ， 一 般 而 言 ， 相 平面 上 的 一 条 闭合 轨道 的 作用 量 按 定 
义 即 为 轨迹 所 包围 的 面积 ， 正 负 取 决 于 定向 。 事 实 上 ， 这 正 是 Henrard 等 人 的 半数 值 摄 动 
方法 的 基本 思想 ， 即 通过 经 典 的 摄 动 理 论 导出 系统 关键 量 (如 作用 量 、 本 征 频率 、 拟 积分 等 ) 
的 分 析 表 达 式 ， 进 而 用 轨道 积分 等 数值 方式 对 这 些 分 析 表 达 式 进行 数值 估计 。 由 式 (47) 得 
到 的 五 及 I， 可 以 分 析 或 数值 地 反 解 出 : 


=Q (L, h, 01) 
ny 158 


=} 


; (48) 
qo = Q> (Io, 05 + p(h, Io, 61)) 
= P (Io, 62+ p(h, I>, 61)) 
其 中 ， 6 
` /AQ, OP; | 
Ld.) do, . 49 
pl 1, 42, 1) J (Gr 86, 80, OL 1 ( ) 
不 难 验证 式 (48) 是 一 个 正则 变换 ， 将 式 (48) 代入 式 (46) 得 : 
K = Ko (Áh, I) +eF (L. To, 0i, 62) š (50) 
利用 可 积 近 似 Ko 对 系统 正规 化 ， 记 系统 无 摄 频率 为 : 
| ÓKo | OK 
Q = 2 R = ms (51) 
定义 泊 松 括号 为 : 
yo Of dg Of Og 
Uis Z 00,01, 01,00; ` (62) 
则 显 形式 正则 变换 ( 记 生成 函数 为 W) 满足 以 下 方程 s 
P (Io Tay 01) = z- f F (I Tay 815 82) 0a (53) 
一 {Ho; W} = = a + a = = F(h, I5, 01, 62) 一 F(h, Io, 01) . (54) 
系统 在 W 定义 的 变换 下 正规 化 到 小 参数 的 一 阶 量 ， 此 时 系统 的 拟 积 分 五 定义 为 : 
L = In —e{;W} +O (e°). (55) 


“pao ea 下 可 以 展开 为 0, 传 里 时 级 数 形式 时 ，Henrard 等 人 给 出 拟 积分 五 基于 式 (55) 的 
数值 计算 方法 ， 有 具体 细节 这 里 不 再 叙述 。 


(C hinay iv 4⁄FHRTI| 
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Henrard 等 人 用 以 上 半数 值 摄 动 方法 研究 了 日 木 系 的 3/1 平 运动 共振 ， 并 研究 了 与 
Wisdom 等 人 摄 动 方法 的 差异 。 对 比 前 面 介绍 的 两 种 方法 ， 可 以 看 到 Wisdom 等 人 利用 3⁄1 
#k 振 的 一 个 单 摆 近 似 式 (41) 直接 得 到 了 拟 积分 1; 而 在 半数 值 方法 中 ， 首 先 由 可 积 近似 Ho 
得 到 作用 量 鳌 ， 进 而 通过 正规 化 将 五 修正 到 摄 动 的 一 阶 量 ， 即 式 (55)。 关 于 两 种 方法 的 对 
比 还 需 提 到 的 两 点 是 : (1) 为 了 得 到 单 摆 近 似 进而 求 角 作 用 量 ， 经 典 摄 动 方法 无 一 例外 地 
在 平衡 点 处 对 无 摄 项 展开 并 截断 至 平方 项 (ECL i rH H H) 导出 HI! 正 是 这 一 处 理 方法 )， 
Henrad 等 人 指出 这 一 做 法 使 得 截断 系统 的 拓扑 结构 在 远离 平衡 点 处 ， 特 别 是 在 分 界线 附近 
可 能 严重 地 偏离 原 系统 ， 而 通过 数值 方法 计算 给 定 轨道 的 角 作 用 量 不 需要 做 此 截断 ，(2) 式 
(51) 即 系 统 无 摄 频率 同样 可 以 用 数值 方法 进行 计算 ， 则 进而 可 以 得 到 O, 与 O, 成 整数 比 ， 
即 次 级 共振 的 位 置 。 

下 面 我 们 简单 讨论 以 上 介绍 的 长 期 摄 动 理论 在 1/1 共振 中 的 应 用 。1/1 共振 中 与 式 (41) 
的 HI 对 应 的 是 式 (39) 的 HI, ipm ELO 节 中 提 到 的 那样 ，Hi1 的 角 作用 量 并 不 能 简单 地 
写 出 ， 故 应 用 绝热 摄 动 理论 ( 即 式 (45)) 去 显 式 地 计算 拟 积分 十 分 困难 ; 需要 指出 的 是 ， 若 
考虑 1/1 共振 的 局 部 性 态 ， 则 可 将 Hi 在 共振 角 的 平衡 位 置 处 展开 ， 此 时 哈密 顿 函数 可 再 
次 简化 为 谐振 子 或 单 摆 ， 进 而 可 应 用 经 典 的 摄 动 理论 ”。 另 一 方面 ， 依 照 式 (44) 沿 着 1/1 
t 振 可 积 近 似 决定 的 菜 一 轨道 对 哈密 顿 函数 或 哈密 顿 正则 方程 进行 积分 平均 是 可 行 的 ， 
若 考虑 系统 平衡 点 处 的 长 期 演化 ， 则 平均 化 退化 为 将 平衡 点 位 置 直接 代入 原 哈密 顿 函 数 ， 这 
一 做 法 在 1/1 共振 中 已 有 许多 应 用 ”"””。 从 角 作 用 量 的 计算 来 看 ，Henrard 等 人 的 半数 值 
摄 动 方法 似乎 十 分 适合 1/1 共振 ， 然 而 ， 该 方法 基于 经 典 摄 动 理论 要 求 摄 动 函数 能 够 展开 为 
1/1 共振 角 的 傅 里 叶 级 数 形式 ， 而 这 一 要 求 对 于 振幅 较 大 的 1/1 共振 是 十 分 困难 的 (I TA, 
区 [3 引 节 )， 如 何 将 Henrard 等 人 的 半数 值 摄 动 方法 应 用 在 1/1 共振 上 是 一 个 值得 研究 的 问 


题 。 


本 节 主 要 介绍 了 1/1 共振 在 无 其 他 摄 动力 下 的 长 期 摄 动 理论 。 当 存在 其 他 摄 动 ， 例 如 其 
RRS JRT, KEE, Poynting-Robertson (P-R) 效应 等 辐射 效应 ™、 行 
星 盘 的 摄 动 或 系统 参数 的 缓慢 变化 时 ”™， 系 统 1/1 共振 态 的 稳定 性 、 稳 定 范围 及 长 期 演 
化 是 一 个 重要 的 研究 方向 。 

2.2 Hill 三 体 问题 
2.2.1 Hill 方程 

记 太 阳 、 月 球 和 地 球 质量 分 别 为 mo, mi, ma; MWA m/m» z 1/80 < 1, u! = mə/(mo+ 
m) z 3 x 10 < 1， 故 月 球 在 地 球 、 太 阳 的 引力 作用 下 的 轨道 运动 可 近似 为 限制 性 三 体 问 
题 。 为 建立 月 球 运 动 的 中 间 轨 道 ，Hill 在 限制 性 三 体 问 题 的 框架 下 进一步 考虑 了 以 下 近似 : 
(1) 月 球 在 黄道 面 内 运动 ， 且 地 球 和 太阳 的 开 普 勒 运 动 为 圆 轨 道 ， 则 系统 简化 为 平面 圆 型 限 
制 性 三 体 问 题 ，(2) 记 月 球 和 太阳 与 地 球 的 距离 分 别 为 x A Ry, W r/R = 2.5 x 107%, WW 
ATLA r/R, 和 jw 为 小 量 对 摄 动 函 数 曙 级 数 展开 ， 并 省 略 所 有 与 7/Ri 相关 的 项 。 考 虑 以 上 近 
似 ， 在 原点 在 地 球 的 日 地 会 合 坐标 系 下 ， 使 用 直角 坐标 ， 系 统 的 运动 方程 退化 为 简单 的 代数 


f^ ° AXE É LAN. ¿(= H oT] 
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微分 方程 中， | 
.. . H | 
X-3y=- T X +3X 
m ’ (56) 
Y +2X = -Rr 


KF, (X, Y) 是 月 球 的 地 心 坐 标 ，R? = X? 十 Y2。 在 经 


文献 中 ， 系 统 (56) 被 称 为 Hill 问 


题 ， 因 式 (56) 是 在 一 定 假设 下 从 限制 性 问题 推导 而 来 ， 故 一 段 时 间 以 来 人 们 普遍 认为 Hill 
问题 是 限制 性 问题 的 一 个 子 集 ， 然 而 Hénon 和 Petit’ 首先 指出 ， 当 
十 
Mı Mm, m&m < cn «1, (57) 
0 


并 且 mi 与 mo 的 相对 昌 


Ea. 
里 


] 于 一 个 天 体 的 质量 远大 于 其 他 两 者 的 情形 ( 即 式 (57))， 
质量 远 小 于 其 他 两 者 的 情形 ， 即 


Mı «m», 


pER 


E 离 在 (w) ERIT, mi, ma 相对 i 
(56) 有 相同 的 形式 ， 由 此 ，Hill 问题 实际 上 与 限制 性 问题 具有 可 


Mı < mo ù 


运动 的 微分 方程 恰好 与 经 典 的 式 
日 同 的 “地 位 ”: Hb 问题 适 
由 性 问题 适用 于 一 个 天 体 的 


而 限 


(58) 


H 5j 


般 情 形 ( 式 (57) 成 立时 ) Hill 方程 的 主要 思路 是 : 在 


个 “平均 ”转动 多 


WARA F, 3⁄4 


M1, ma 的 运动 解 耦 为 mi, ma 质心 相对 于 mo 的 运动 以 及 mi, ma 的 相对 运动 。 按 “平均 ” 
转动 坐标 系 为 圆 轨道 或 椭圆 轨道 之 分 ，Hill 方程 也 有 圆 型 及 椭圆 型 之 分 ， 下 面 给 出 平面 椭圆 
型 Hil JFE. ERF, m, ma 相对 于 中 心 天 体 mo 的 位 置 矢量 ri, rs 满足 微分 
方程 : 
" d 
Pı = — (1 + em) ne 十 EH2 8 
Ir; | Ira ld 
(59) 
d 
To = — (1 + €p2) A + Ep — uc 
[rə] Ir| Id 
其 中 通过 对 单位 归 一 化 ， 使 得 G —1, m=1, 3FH. 
He, eo. hao BS. (60) 
mo m; + mo 
引入 mi, m» 的 质心 : 
maipi 十 Mor 
R= = = M171 + M272 , (61) 
Mı + M2 
不 难得 到 以 下 关系 : 
rı = R- pod, T; = R+ pmd (62) 
将 式 (62) 代入 式 (59)， 不 难得 到 : 
| Hı a H2 2) 
Iri] [r2] (63) 
= T1 T» d 
Ir? dr Ja}? 
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考虑 m,m 相互 距离 为 小 量 ， 即 令 


O(R)=1, O(d) < 1 . (64) 


此 时 ， 不 难得 到 以 下 展开 式 : 


[r1] T ee 4 o (d?) 
R R5 (65) 
一 3 1 ui Rd 2 
Ira| = ps ? m + O (d?) 
Jü, R.d 为 矢量 的 点 乘 ， 将 式 (62), (65) 代入 系统 (63)， 可 得 
R= (14) (- o0) 
i (66) 
“ d d (RAR 
a=- m^ m ^9) 
^ 
d= r => d=el/3r, (67) 
将 式 (67) 代入 系统 (66)， 可 得 
R 
— TOEN] " 
T T (R-r)R : 
fone Ir]? R3 F3 RS | O (e?) 


在 给 出 进一步 的 推导 之 前 ， 我 们 对 式 (64) 和 式 (67) 做 简单 解释 : 式 (64) 的 近似 表明 系统 
(66) 适用 于 描述 mi, m 相互 交会 时 系统 的 运动 ， 并 且 此 时 系统 的 质心 运动 从 完整 运动 方 
Ke PAA; Fst (67) 给 出 量 级 估计 ， 因 此 系统 (68) 截断 至 小 量 d, e 的 零 阶 时 保留 了 最 
多 的 项 ， 且 此 时 系统 不 再 含 任何 参数 ， 从 物理 上 看 ，d = cU? IE Hill 半径 的 定义 ， 此 时 
ma, ma 的 相互 引力 作用 与 mo 的 中 心 引 力 在 同一 量 级 ， 故 此 时 运动 方程 截断 至 小 参数 的 零 
阶 项 时 ( 即 系统 (68)) 保留 了 两 种 相互 作用 的 零 阶 项 。 下 面 继续 推导 椭圆 型 Hill 方程 ， 为 简 
化 记号 ， 将 二 维 实 矢量 R, r 视 为 复 变量 RR, r e€ C， 不 难 验 证 ， 矢 量 的 点 积 Rr 在 复数 运 


算 下 可 表示 为 : 
= ar 
R.r=Re(R P) ， (69) 
其 中 ，Re 代表 求实 部 。 为 了 在 随 质 心 R 一 同 运动 的 脉动 会 合 坐 标 系 下 考虑 相对 运动 ， 引 入 
复 变量 z, 
r = Rz (70) 


由 式 (70)， 不 难 推 得 z 的 含义 为 脉动 会 合 坐标 系 下 m 和 m 的 相对 运动 矢量 的 复数 形式 。 
将 式 (70) RAR (68) 的 第 二 式 ， 结 合式 (68) 的 第 一 式 及 式 (69)， 不 难 验证 
d?z 1dRdz E 3 


+ 2 = Fz Rez. 71 
d? R dt dt |Rz| F? = SH 
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A (71) 为 脉动 会 合 坐 标 系 下 的 相对 运动 方程 。 由 式 (68) 的 第 一 式 ，mi, ma 质心 近似 作 开 普 
勒 运 动 ， 则 有 以 下 关系 : 


. a (1 — e?) ds a (1 — e2) 
= = = . 72 
HOO) ae di l+ecoss’ dt R? (2) 


为 了 进一步 简化 系统 (71)， 做 以 下 时 间 尺 度 变 换 : 将 对 时 间 t 的 求 导 变换 为 对 真 近 点 角 s 的 
求 导 ， 由 式 (72) 不 难 验证 以 下 关系 : 


d d q? esins a (1 — e?) 
=A B A=-2 B= 
dt? ds " ds? ' R 7! R4 ; (73) 
B-1A —2esins 1 dR |. 2esins es 1 1 
 l+ecos ' “Rds l+ecoss ` ' BR) 1+ecoss 


将 式 (73) 代入 式 (71)， 可 得 : 


d2z dz 1 z 
21 | 3R =0. 74 
ds? ' kam (4 ez) vet) 


式 (74) 即 为 椭圆 型 Hill 方程 。 令 e—0,z-—zc--iy, WAR (74) 退化 为 圆 型 Hill 方程 : 


$—2j4 (a? +y?) "2-32 = 0 
ü+2ë+ (a? +) 

不 难 发 现 ， 系 统 (75) 与 经 典 的 Hill 方程 ( 见 (56) st) 有 相同 的 形式 ， 且 不 含 任何 参数 ， 由 于 
不 含 任何 参数 ， 当 mi 与 ma 很 小 时 (例如 接近 或 小 于 计算 机 的 机 器 精度 时 )， 系 统 (74) 或 
(75) 比 完整 的 三 体 运动 方程 (例如 系统 (59)) 更 适合 数值 积分 ” 。 
2.22 平面 圆 型 Hill 方程 的 解 

系统 (75) 存在 两 个 平衡 点 (3-73 ,0)， 分 别 对 应 于 圆 型 限制 性 三 体 问 题 中 的 Ly D» 
点 ; 系统 (75) 并 不 存在 与 Ls 点 对 应 的 平衡 点 ， 原 因 在 于 推导 过 程 中 式 (64) 所 作 的 假设 ， 以 
及 在 推导 中 去 掉 了 所 有 与 视差 项 (d/ R) 相关 的 项 ， 因 此 在 以 上 假设 与 近似 下 ，mo 以 及 Ls 
点 的 位 置 在 -co 处 ， 自 然 无 法 出 现在 系统 (75) 中 。 类 似 于 圆 型 限制 性 三 体 问 题 ， 贺 型 Hill 
方程 ( 见 式 (75)) 存在 形式 如 下 的 运动 积分 h: 


(75) 


h= 5 (Pg) - (a ty? os. (76) 


系统 (75) 在 一 定 近似 下 可 以 给 出 分 析 解 ， 下 面 介 绍 相关 内 容 。Dermott 和 Murray” 对 圆 型 
限制 性 三 体 问 题 中 的 马蹄 型 轨道 进行 了 系统 的 数值 研究 ， 他 们 将 小 天 体 ma 放 在 与 ma 的 初 
始 相 位 相差 180° 的 位 置 Bn 固定 me 的 轨道 半 长 径 ， 在 不 同 的 初始 半 长 径 a, 下 积分 mi 的 
轨道 ， 数 值 积分 发 现 当 ar 全 ai — as RETER, m 运动 在 几乎 对 称 的 马蹄 型 轨道 上 。 下 
面 从 Hill 方 程 角度 研究 这 一 情形 ， 进 而 从 系统 (75) 导出 一 个 具有 类 似 特 征 的 解 。 从 物理 直观 
LG, hs 0 意味 着 初始 zo 0; m 与 mo 的 初始 轨道 相位 差 为 180"， 从 在 Hill 方程 的 近 
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似 式 (64) 下 意味 着 初始 yo © oo; 而 一 条 几乎 对 称 的 马蹄 型 轨迹 在 Hill 坐 标 (z, y) 下 表现 为 
一 条 关于 z = 0 对 称 的 U 型 轨迹 ， 若 zo > 0， 则 + Pi, Ty 的 符号 与 z 符号 相反 ， 
可 以 猜测 yx 一 xz。 由 以 上 分 析 ， 在 系统 (75) 第 一 式 中 令 (a? +y?) — 0, 2 — 0, W 


2 -3 (77) 
将 式 (77) 代入 系统 (75) 的 第 二 式 ， 得 : 
2 
之 一 CE (78) 
将 式 (77) 与 式 (78) 相 除 ， 分 离 变 量 后 积分 得 到 : 
8 
y= 3(ci z2) : (79) 
式 (79) 给 出 的 正好 是 一 条 U 形 轨迹 ， 且 运动 积分 c 即 为 y 3 oo IE) “MBSR” he 
数值 积分 发 现 ”， 当 c < 0.7 时 ， 式 (79) 给 出 了 原 系统 ( 见 式 (75)) 数值 积分 轨迹 非常 精确 


的 近似 。 
当 m, ma 距离 充分 大 以 至 于 它们 的 相互 引力 可 以 忽略 时 ， 即 (a2 +y?) 3? 一 0， 式 
(75) 简化 为 : 


ë — 2 — 3r = 0 
. (80) 
ë+ 2 =0 
系统 (80) 是 一 个 常 系数 齐 次 微分 方程 组 ， 其 通 解 为 : 
多 = C, + Co cos (t — tı) 
y= ced (t = t*) E 2C» sin (t = tı) 


Rh, Ci, Co, €, ti 是 运动 积分 ， 事 实 上 ， 这 些 运动 积分 有 明确 的 物理 意义 ， 因 为 系统 (80) 
忽略 my, ma 的 相互 引力 ， 此 时 mi, mo 作 固定 的 开 普 勒 运动 ， 故 式 (81) 为 Hill 方程 的 无 摄 
解 ， 此 时 可 以 建立 Ci, Co, t*, ti 与 mi, mo 轨道 根 数 的 函数 关系 ， 相 关 推 导 可 见 区 ZZ 节 。 

X (81) 对 应 着 上 一 oo 时 系统 的 解 ， 由 式 (81) BZR yx t Mest (81) 也 可 视 为 y 一 oo 
时 系统 的 解 。 一 个 自然 的 想法 是 : 能 否 找到 一 组 以 1/t 或 1/y 为 小 参数 的 容 级 数 解 ， 使 得 当 
t — oo (或 y — oo) 时 ， 该 寡 级 数 解 渐 近 逼近 于 解 (81)。Spirig 和 Waldvogel” 以 及 Hénon 
和 Petit 分别 给 出 了 一 族 四 参数 渐 近 逼近 解 ， 以 1/t 为 小 参数 的 解 的 形式 如 下 : 


oo j oo 3 
= 一 k = 1— k E 
r= 5 ti > ajk (t), y= 5 i > bg (t). ilg "E (82) 


j=0 j=0 


Ah. ajn(t), b; (t) 是 时 间 上 二 的 三 角 函 数 多 项 式 (包括 常数 )。 可 以 看 到 ， 式 (82) JPR EH 
见 的 震级 数 解 ， 而 是 包含 了 时 间 t 的 对 数 项 ， 这 使 得 上 一 co 时 ， 式 (82) 实际 上 并 没有 收敛 


— 
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T Hill 方程 无 摄 解 ( 即 式 (81)). Hénon 和 Petit ™ AH T A 1/y 为 小 参数 的 渐 近 逼近 解 ， 该 
解 同 样 包含 y 的 对 数 项 。 前 文 已 经 提 到 ，Hill 方程 的 无 摄 解 可 以 看 作 系 统 的 外 部 近似 ， 即 
当 mi 与 m2 的 相对 距离 与 m1, ma 质心 到 mo 的 相对 距离 在 同一 量 级 时 的 系统 运动 的 近似 ， 
而 当 m 与 ma 近 距 离 交 会 时 ， 完 整 Hill 方程 的 数值 积分 可 看 作 原 系统 运动 的 高 精度 的 内 
部 近似 ， 而 式 (82) 给 出 的 渐 近 逼近 解 恰 好 是 内 部 近似 与 外 部 近似 的 E, Hill 77 
程 被 用 来 处 理 一 类 称 之 为 “卫星 遭遇 ”(satellite encounter) 的 问题 ”。 这 一 类 问题 的 特点 
是 : m 与 m 初始 距离 较 远 ， 由 于 半 长 径 不 同 而 在 空间 中 人 靠近， 在 有 限 的 时 间 内 处 于 近 距 
离 的 相互 作用 阶段 ， 然 后 在 相互 引力 作用 下 再 次 互相 远离 。 以 圆 型 Hill 方程 为 例 ， 应 用 式 
(82) 研究 “卫星 遭遇 ”问题 的 具体 方法 是 : 式 (82) 中 的 四 个 运动 积分 对 初始 状态 空间 参数 
Wo m 与 m» 的 初始 靠近 阶段 由 式 (82) 给 出 ; 当 mi 与 mo 的 距离 小 于 某 个 六 值 时 OZBE 
取决 于 式 (82) 的 精度 ， 进 而 取决 于 式 (82) 的 阶 数 以 及 设 定 的 误差 闵 值 )， 系 统 的 运动 由 完 
整 Hil 方 程 数值 积分 给 出 ， 在 相互 引力 作用 下 当 m 与 ms 的 相互 距离 大 于 给 定 阐 值 时 ， 系 
统 的 运动 再 次 由 渐 近 解 给 出 。 


= 2.2.3 圆 型 Hill 方程 的 正则 ie 
< 在 哈密 顿 系 统 中 我 们 看 到 ， 通 过 定义 合适 的 正则 变量 ( 即 定义 1/1 共振 角 等 )， 系 统 不 同 
N B occa. 进而 可 以 使 用 正规 化 或 平均 化 消除 系统 的 快 变量 。 


wies, AA Hill 方程 的 无 摄 解 可 以 与 ma, ma 的 开 普 勒 根 数 建立 联系 ， 而 开 普 勒 根 数 可 
进一步 转换 为 Delaunay 根 数 等 正则 变量 ， 这 启发 我 们 通过 Hill 方程 的 无 摄 解 引入 mi, ma 
的 开 普 勒 根 数 ， 进 而 对 ma 和 m 的 相对 运动 定义 一 组 合适 的 正则 变量 ， 阁 这 组 变量 同样 具 
有 等 级 式 的 特征 (事实 上 的 确 如 此 )， 则 可 类 比 于 哈密 顿 系 统 对 Hill 方程 或 Hill 方程 的 能 量 
积分 进行 平均 化 消除 快 变量 ， 进 而 研究 平均 化 系统 的 动力 学 性 质 。Namouni” 在 工作 中 首先 
给 出 了 Hill 方程 的 正则 形式 ， 下 面 对 这 一 形式 作 简 单 介 绍 。 考 虑 三 维 情形 的 圆 型 Bil 方程 ， 
即 在 系统 (75) 中 加 入 Z 方向 的 运动 ， 可 得 : 


£ — 2g + (z? + 1° 2) r- 3z = 0 
j + 22 + (2? y? 2) y=0 . (83) 


E T 


Z+ zc(z + 2 + 22) 
空间 圆 型 Hill 方程 ( 见 系统 (83)) 的 无 摄 解 为 : 


— di 
v= SO G= t") —2Cysin(t—th) , (84) 
ee 


HH, S= (Ci, Cz, C3, t*, t, to} 是 运动 积分 ， 下 面 我 们 建立 S 与 mi, M2 mr: AR 
的 函数 关系 。 由 式 (62), (67) 及 (70) 得 : 


rı = R(1 — me"? z) ] r2 = R(1 + me"? z) ; (85) 
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此 时 z = (z, y, z) H Hill 方程 无 摄 解 (3X (84)) 给 出 ， 则 
Rn —1 Sel? C,-— josef? C» COS (t — t1) + uE” a (t — t* )+ 
912el/3 C^» sin (t = tı) 一 /2E 1⁄3 C3 sin (t = 12) š (86) 
R r 有 类 似 的 表示 式 ， 这 里 不 再 列 出 。 式 (86) 中 Rr 代表 m 在 由 质心 运动 (B R) 所 
决定 的 脉动 会 合 坐 标 系 下 的 运动 。 事 实 上 ， 玉 -7 可 视 作 m 相对 于 “本 轮 ”( 即 R) 所 作 的 
“ 均 轮 ”运动 (epicyclic motion)。 由 于 所 考虑 的 是 圆 型 Hill 方程 ， 此 时 R 对 应 一 个 开 普 勒 
轨道 ， 记 该 圆 轨道 的 轨道 半 长 径 为 w%， 轨 道 角 速度 为 必 ， 因 mo = 1,mi,m2 < mo, A 
n! = 2/7, W m 的 轨道 半 长 径 为 a;， 则 由 式 (64) 得 : 
O(a; — a’) = O(e;) = O(¿;) = € cosi; 1, sini; ei; , (87) 


其 中 ，j = 1 2, ej, i; 分 别 为 m; AE m ty ER UTE TAA. FILA PEAR RC UB eE e 
m 在 本 轮 及 下 的 均 轮 运动 表达 式 ， 在 一 个 以 mw 为 角速度 的 匀速 转动 坐标 系 下 ， 考 虑 到 式 
(87), m, 的 位 置 矢 量 Pi 为 : 


Tzu 


rı cos (fı + à — n't) 
Ti = R. (n't) Tı = Ti sin (fi + Wy = n't) n (88) 
rı cos (fy + w1) 


RP, fi, w, Qi, Qi = wi + Q) DAE mi 轨道 的 真 近 点 角 、 近 点 角 距 、 升 交点 经 度 和 近 点 
角 经 度 考虑 到 O(e;) = O(i;) = e， 对 式 (88) 右 端 进行 椭圆 展开 并 截断 至 一 阶 项 ， 得 : 


>. 


El 


1 
(cos (Mi—n't+ar) +51 (cos (2M, —n't+@1)—3 cos Oo) 
1 
Tı =a, (sin (M, —n't+ar)+5e1 (sin (2M, —n't-à1)4-3 sin (n't-a))) +O (&) : (89) 


1 
21 (e (sin (2M, cu) —3 sin w1) +sin (Mı +w) 


ER (64) KRF, RIIA w (t-t) — (M +Q) = 0 (e)， 其 中 要 取决 于 空间 惯性 系 
X ANH, 4 t 一 0， 则 式 (89) 可 简化 为 ; 


al 一 alel cos (n't — à) 
Tı = | aa (M, — n't +a) + 2aje;sin(n't-@)] ` (90) 
71a, sin (n't — £21) 
取 本 轮 有 R 圆 轨道 的 轨道 半 长 径 为 单位 长 度 ， 则 有 a =n’ = 1， 此 时 ， 由 式 (86) 与 式 (90) 分 
别 给 出 的 m 的 均 轮 运动 的 差别 应 该 为 一 阶 小 量 ， 即 Rr- r = O (e)， 由 此 解 出 : 
| mer Ci =1-a,, yigg Co = a161 , we Cs = —iia1 


3 EM (91) 
50 -a)t = —(21 , tə = Qı 


ti =w, 
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对 ma 有 类 似 的 关系 ， 这 里 不 再 列 出 。 将 式 (91) RER (84)， 即 得 到 Hill 坐标 在 m, 轨道 根 
数 下 的 表达 式 ， 类 似 地 可 得 到 m» 轨道 根 数 表示 的 Hill 坐标 ， 这 里 不 再 列 出 ， 结 合 mi, me 
的 表达 式 ， 不 难得 到 : 


Z = hT HUT & [az 一 ai 十 alel cos (t(—@1) — azez cos (t—W)| 


ETT 


1 [3 : MEM — 
Y = IV Ho © as 5 (a9 —a4) (t—t*) 4-2a3es sin (t—W2) —2a ,e; sin(t—@,)| , (92) 


1 
Z = Zt Hz = ae [io sin (n't— (25) — i, sin (n't— £21)] 


下 面 我 们 由 mi, mo 的 轨道 根 数 定义 mi, m 相对 运动 的 轨道 根 数 。 不 难 验证 ， 由 极 坐 标 式 
的 正则 变量 (I, 0) 到 直角 式 坐 标 (p = V —21sin0, q = /—2Icos0) 是 一 个 正则 变换 ， 故 从 
无 奇 点 形式 的 Delaunay 变量 ( 见 式 (8)) 可 以 定义 如 下 正则 变换 : 


(L, àn Gi, Gi, Hh, 2) > (Ly, Ai, hi, ki, Pi, d), (93) 
pus hy =V Aen sin à = CE1 Sind, 5 ky =y A SN COS Wy = €] COS (1 P 
(94) 
Pı = —2H, sin f2, = i, sin f, , Gi = V —2R, cos (21 ~ i; cos (2, . 
类 似 地 ， 可 对 ma 定义 如 上 变量 。 由 式 (94), E X mi, m2 相对 运动 的 正则 根 数 如 下 : 
hyp =hg—hy, k.—k;—kÀi, p.—po-m, q=q Q. (95) 
不 难 验 证 ， 使 用 式 (95) 定义 的 相对 运动 根 数 ， 式 (92) 可 简化 为 : 
1 " 
= = |a, — er cos (t — à.)] 
1 
U = um [N + 2e, sin (t — Q,)] ， (96) 
= An sin (t — 12.) 
= Hh, a, = az-a, À, = Aə2— Aro IN (96) 及 其 微分 式 给 出 了 Hill 化 标 及 速度 矢量 (x, y, z, £, 
一 9, 2) 与 相对 运动 根 数 (ar, Ar, hi, krs Drs q.) 的 转换 关系 ， 这 一 几何 关系 可 类 比 于 二 体 运动 
FI 中 的 位 置 速度 与 轨道 根 数 的 转换 关系 。 由 圆 型 Hill 方程 ( 见 式 (83)) 不 难 推导 相对 运动 根 数 


满足 如 下 的 微分 方程 : 


Ow ; Ow . Ow 
Ar — — , hy = —— š Dr = — 

al, Ok, ðq (97) 
i 1 n Ov n Qv ] Ow ^" 
r= une , — ui r = 

2 Oar ðh, s: opr 
其 中 ， ls E 和 r/2， 
p = —-2/3 (x? +y? + z2) 
(98) 


—1/2 
mox [c — e, cos (t — &,))? + (X, + 2e, sin (t — &,))? + isin? (t — 2.) 


(O 


k 
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在 相对 运动 根 数 下 ， 了 ill 方程 的 运动 积分 变 为 : 


1 3 
H,=2h= (2 + i2 a) Ly, (99) 


可 以 看 到 ， 式 (97) 表现 出 与 哈密 顿 系统 类 似 的 正则 形式 (但 需 注 意 系统 (97) 并 不 是 一 个 哈 
密 顿 系统 )， 故 不 妨 称 式 (97) 为 圆 型 Hill 方程 的 正则 形式 。 在 相对 运动 根 数 的 基础 上 ， 可 定 
义 如 下 的 正则 变换 ， 将 变量 重新 变 为 直角 型 坐标 : 


h.. k J, = —, Wr = Wr — 2, 
Ta MTS 2 
4 2 2 ? 
Pre) (ki a (100) 
r 2 > der 


dJ, A dw, +dK, A dQ, = dh, ^ dk, + dp, ^ dq, . 


2.2.4 Namouni 的 摄 动 理论 
从 Hill 坐标 与 相对 根 数 的 变换 关系 ( 见 式 (96)) 可 得 ， 系 统 的 运动 可 视 为 自由 度 (ar, Ar) 
的 运动 与 一 个 周期 摆动 的 复合 ， 其 中 周期 摆动 在 zy 平面 内 的 振幅 由 相对 偏心 率 er 决定 ， 
周期 对 应 于 系统 的 平 运动 时 标 ，z 方向 短 周期 摆动 的 振幅 由 相对 轨道 倾角 i HR; Al 入 
正 是 哈密 顿 系统 中 的 1/1 共振 角 ， 故 一 个 自然 的 假设 是 入 的 时 标 远 大 于 系统 平 运 动 t 的 
时 标 ， 此 时 若 a, > e:， 则 系统 的 运动 表现 为 (a, 入) 决定 的 主导 运动 和 短 周期 振动 的 复 
合 ，Namouni 的 研究 正 是 基于 这 一 自由 度 的 等 级 式 结构 。 事 实 上 ， 关 于 自由 度 的 等 级 式 结 
构 还 需 说 明 er, dr f£ (ap, À n 的 时 标 下 变化 十 分 缓慢 (绝热 变化 )。 在 平面 圆 型 Hill 问题 下 ， 
Hénon 和 Petit 用 Kruskal ”的 摄 动 方 法 证 明了 式 (82) 的 两 个 参数 (碰撞 参数 和 约 化 偏心 
A) 在 共 轨 运动 的 时 间 尺 度 下 是 绝热 不 变量 ， 而 碰撞 参数 对 应 上 一 co 时 a. 的 值 ， 约 化 偏心 
率 正 是 eps MAF ir 的 绝热 变化 尚 无 分 析 上 的 证 明 ，Namouni 通过 数值 积分 计算 发 现 ， 对 
于 系统 (97) 的 规则 运动 ， 关 于 i, 的 假设 是 合理 的 。 
根据 以 上 分 析 ， 可 以 对 Hill 方程 的 正则 形式 在 一 个 平 运动 周期 ( 即 (0,2z)) 内 积分 平均 ， 
进而 在 平均 化 系统 中 将 (er, ir, w;:， 人 2;)“ 冻 结 ”， 此 时 平均 化 系统 退化 为 自由 度 可 积 系 统 ， 
其 轨迹 由 平均 化 的 运动 积分 给 出 : 


UL) = K, —-a + (W) , (101) 


à, l, 2 MAE. š 
"i. COS u ui —+sinu} + —1 | sin? (ut+w,) 
27€, er e2 
0 


HH, u-t-G. AW) USA wo H (101) WA, K, 是 平均 化 系统 的 一 个 运动 积 4 
“1, 一 co 时 (Bl mi, ma 的 相互 距离 与 质心 距 R fE[R] RARE), (0) 一 0， 此 时 碰撞 参数 
ao 可 定义 为 : 


1 


du, (102) 


2 8 


Qo = 3 (K, = H,) Y (103) 


EH 
H 


H 


F! 
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H3X (101) 可 知 ， 系 统 主导 运动 的 轨迹 满足 如 下 方程 : 


8 G l 2K. 

2 2 r T T 

a lr =a S 》 , 1 ; Wr , 
u ( ) 9 3er (s er ( e2 ) ) 


juni 


2m -i (104) 
S (a, 1, R, w) = = f [à cos + 4(T+ sinu)? + (R 1) sin? (to)] du 
0 
注意 到 方程 左 端 和 右 端 均 显 含 a,，Namouni 进而 考虑 a, = 0 时 共振 角 刀 的 取 值 0， 则 有 : 
m. - l 2K, 
gaoer = S (° "* d vr) i (105) 
Bl a2 (Ip) > 0， 故 在 a, == 0 BF, 1, 的 取 值 必须 使 得 
l, 2K, lo 2K, 
S (° e." e P w) < S (o, e e 5 3 x (106) 
- Namouni 定义 的 函数 5 很 好 地 描述 了 平面 及 空间 情形 下 系统 m, ma 相对 运动 的 可 能 
= 类 型 ， 这 一 方法 可 与 哈密 顿 系统 的 一 自由 度 “ 冻 结 ” 相 图 ( 见 图 团 进行 比较 。 为 研究 系统 的 
N 长 期 演化 ，Namouni 进一步 对 正则 形式 的 Hill 方程 积分 平均 ， 利 用 平均 化 的 摄 动 方程 研究 
C^ 了 不 同 共振 态 下 轨道 进 动 的 频率 及 长 期 共振 等 性 质 ， 其 思想 与 哈密 顿 系统 中 的 平均 化 方法 


没有 本 质 上 差异 ， 故 不 再 叙述 。 

由 式 (106) 即 可 确定 1, 的 运动 “趋势 ”， 图 加 给 出 了 两 组 不 同 参数 下 函数 S 的 取 值 (5 关 
于 7 是 侦 函数 ， 故 只 画 出 7 > 0 情形 )。 图 Ba) 为 平面 情形 i = 0， 此 时 按 定义 玉 =1,，17=1 
是 函数 S 的 奇 点 ， 当 lo 取 值 在 图 中 蓝 色 虚线 与 实 线 的 交点 时 ， 则 按 以 上 讨论 有 1 > do 3X 
| < 一 l0， 故 这 一 区 域 对 应 马蹄 型 轨道 ， 类 似 地 ， 红 色 曲 线 与 实 线 右 支 的 交点 也 对 应 着 马蹄 
型 轨道 ， 而 当 lo 取 在 红色 虚线 与 实 线 左 支 的 交点 处 时 ， 由 以 上 讨论 有 一 lo < 1 < /0， 则 这 一 
区 域 对 应 着 有 界 运 动 ， 从 运动 轨迹 上 看 可 类 比 于 拟 卫 星 轨 道 ， 黄色 虚线 代表 碰撞 参数 ao 的 
较 大 取 值 ， 此 时 共振 角 A 没有 进入 共振 ， 故 称 之 为 穿越 轨道 。 图 邯 ) 为 空间 情形 K = 1.5， 
按 定 义 er < 办， 不 同 颜色 的 实 线 对 应 参数 w, 的 取 值 ， 可 以 看 到 ， 相 对 于 图 a)， 绿 色 虚 线 标 
出 了 一 种 新 的 运动 类 型 ， 当 io 取 值 为 绿色 虚线 与 黄色 实 线 的 交点 处 时 ，! 既 可 能 向 大 于 lo X 
向 运动 ， 也 可 能 往 相 反方 向 演化 ， 当 考虑 ej, i, 的 长 期 演化 时 ， 轨 道 可 能 从 有 界 轨道 跃迁 为 


~] 


u 马蹄 型 轨道 (反之 亦 然 )，Namouni 将 这 一 种 空间 情形 下 的 轨道 类 型 称 为 跃迁 轨道 。 


23 ”数值 计算 

对 完整 系统 运动 方程 进行 数值 积分 是 研究 限制 性 或 行星 三 体 问题 最 直接 、 也 是 最 精确 
的 方法 。 为 减 小 数值 计算 负担 ， 基 于 平均 化 的 方法 ， 可 对 原 系统 构建 一 个 保留 长 期 动力 学 
演化 的 映射 模型 ”™“， 对 哈密 顿 系统 构造 映射 模型 的 标准 方法 可 参考 Hadjidemetriou 的 工 
作 呈 。 将 轨道 数值 积分 与 一 些 数值 方法 相 结合 ， 可 对 系统 相 空 间 结构 作 精 确 的 探索 ， 其 中 
常用 的 数值 方法 有 寻找 和 延 拓 周期 轨道 、 庞 加 莱 截 面 、 稳 定 指标 计算 等 。 类 似 于 平衡 点 在 平 
均 化 系统 中 的 重要 性 ， 周 期 轨道 在 系统 全 局 的 动力 学 结构 中 扮演 着 “骨架 ”的 角色 “， 周 
期 轨道 的 分 叉 、 附 近 的 稳定 及 不 稳定 流 形 决 定 了 相 空 间 局 部 及 全 局 的 性 质 。 因 此 ， 数 值 延 
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1.8 r 
—o,=0° 
| — a =30° 
1.4} — o,—60^ 
— @=90° 
--- transition | | 


a) b) 


ik: a) 平面 情形 ， 其 中 蓝 色 虚线 处 在 马蹄 型 轨道 区 域 ， 红 色 虚 线 处 在 马蹄 型 轨道 或 有 界 轨道 ， 黄 色 虚 线 对 应 穿越 
轨道 ，b) er < i, 情形 ， 绿 色 虚 线 对 应 一 条 可 能 的 跃迁 轨道 。 


5 ”不 同 参数 下 5 (a, l, K, wr) 的 图 像 


拓 周 期 轨道 是 研究 系统 相 空 间 结构 的 常用 手段 ”"”。 庞 加 莱 截 面 "将 维 的 相 空间 和 

N — 2 维 相 空 间 联系 起 来 ， 故 可 直观 地 从 低 维 空间 的 几何 上 “观察 ”系统 的 动力 学 结构 。 同 
时 ， 相 空 < 间 不 同 区 域 的 稳定 性 可 以 通过 稳定 性 指标 来 确定 ”， 例 如 Lyapunov 指数 、 
MEGNO 指标 ”、 本 征 频 率 的 漂移 (NAFF) 等 。 


3 动力 学 简介 


Brown™ E 1911 年 指出 : 在 平面 限制 性 三 体 问题 中 ， 当 质量 参数 pp 足够 小 时 ， 
在 这 样 一族 轨 道 ， 该 族 轨道 起 始 于 围绕 点 L, 的 无 穷 小 振动 ( 即 点 L, 的 长 周期 轨道 族 )， 
着 以 中 心 天 体 mo 为 圆心 ， 半 径 为 1 的 圆 扩张 ， 经 过 点 Ls 处 的 同 宿 轨道 后 演化 为 ane 
道 ， 最 终 终结 于 点 Zi 和 点 Lo. AS PAA GE AES LS 的 同 宿 轨道 ( 即 治 时 间 正 向 
或 负 向 均 趋 于 点 La) 而 演化 为 马蹄 型 轨道 ，Brown” 首先 给 出 了 共 轨 运动 极 坐标 的 一 阶 解 ， 
基于 该 解 Brown 发 现 日 木 系 中 共 轨 运动 的 轨迹 可 以 由 零 速 度 曲 线 很 好 地 近似 ， 而 零 速 度 曲 线 
显示 在 点 Ls 处 绕 转 La, Ls 的 出 时 型 轨迹 演化 为 绕 转 Da, La, Ls 的 马蹄 型 轨道 ，Brown 进 一 
步 计算 了 点 Ls 附近 的 线性 化 流 形 ， 并 发 现 Ls 附近 的 线性 双 曲 流 形 正好 对 应 了 在 Ls 附近 的 
星 蚜 型 及 马蹄 型 轨迹 ， 基 于 该 一 阶 近 似 结果 ，Brown 猜测 完整 的 限制 性 问题 中 有 相同 的 结 


论 。 


在 Brown 之 后 ， 许 多 研究 从 数值 和 分 析 上 证 明了 马蹄 型 轨道 族 的 存在 性 。 然 而， 
Brown 的 关于 点 La 长 周期 轨道 族 经 过 点 Ls 的 同 宿 轨 道 演化 为 马蹄 型 轨道 的 猜想 先后 被 
Deprit 和 Henrard” ”的 数值 结果 和 Garfinkel 的 理论 分 析 ” 否定 。Henrard” 通过 数值 计 
算 发 现 : 日 木 系 限制 性 三 体 问题 的 点 La 的 稳定 和 不 稳定 流 形 并 不 相交 ， 即 Brown 猜测 的 点 
Ls 处 的 同 宿 轨道 对 于 日 地 系 并 不 存在 。 根 据 已 有 的 数值 结果 ”，Henrard 猜测 : 日 木 系 中 
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长 周期 轨道 通过 一 族 周期 轨道 与 绕 转 13 圈 的 短 周期 轨道 ( 记 这 条 轨道 为 13 分 又 轨道 ) 相连 ， 
而 这 条 13 分 又 短 周期 轨道 又 通过 一 族 周 期 轨道 与 14 分 又 短 周 期 轨道 相连 ， 如 此 继续 ， 即 天 
分 叉 短 周期 轨道 与 (k+-1) 分 又 短 周 期 轨道 相连 (k > 13)。 因 此 长 周期 轨道 族 从 三 角 平 动 点 开 
始 ， 依 次 经 过 13 分 又 、14 分 又 、15 分 又 等 短 周 期 轨道 。 由 于 短 周 期 轨道 族 最 终 终结 于 一 条 
绕 转 Ls 的 周期 轨道 ， 因 此 长 周期 轨道 最 终 终 结 于 一 条 绕 转 无 穷 圈 的 短 周期 轨道 (可 认为 是 
绕 转 点 L. 的 一 条 同 宿 轨 道 ， 该 同 宿 轨道 渐 近 于 一 条 绕 转 L, 的 周期 轨道 )。Garfinkel 通过 
1:1 共振 极 坐标 的 二 阶 解 验证 了 Henrard 的 猜想 ， 由 此 从 分 析 上 否定 了 Brown 猜想 。 

人 们 在 太阳 系 内 发 现 的 第 一 个 马蹄 型 共 轨 系 统 是 土星 的 两 个 卫星 Janus 和 Epimetheus, 
关于 Janus 和 Epimetheus 的 观测 历史 可 见于 参考 文献 sl。 根据 观测 数据 可 知 ， 两 个 卫星 
的 质量 比 大 约 是 1:3 至 1:5，Janus 和 Epimetheus 的 半径 之 和 大 约 是 它们 轨道 半径 之 差 
的 3 倍 。 由 于 相互 引力 作用 ， 两 个 卫星 避免 了 碰撞 并 分 别 运 行 在 两 条 马蹄 型 轨迹 上 。 限制 性 
模型 并 不 适用 于 Janus 和 Epimetheus 组 成 的 共 轨 系统 ， 由 此 ，Dermott 和 Murray ™ 38 
过 研究 两 个 卫星 轨道 相互 引力 对 半 长 径 的 摄 动 ， 构 造 了 一 个 适用 于 行星 三 体 问题 的 共 轨 模 
型 。Yoder 5& A" 使 用 定义 在 第 一 主 天 体 上 的 极 坐标 ， 对 非 限 制 性 问题 构建 了 一 个 1/1 理想 
振 模 型 ， 并 根据 该 模型 指出 ，Janus 和 Epimetheus 组 成 的 共 轨 系统 目前 十 分 稳定 。 Spirig 
和 Waldvogel"" 使 用 Hill 三 体 问题 模型 准确 预测 了 Janus 和 Epimetheus 的 交换 轨道 ， 并 对 
系统 的 稳定 性 得 到 相似 的 结论 。 

前 面 提 到 ，Hil 方程 可 用 来 研究 一 类 称 之 为 “卫星 遭遇 ”的 问题 。Hénon 和 Petit” 对 
“卫星 遭遇 ”问题 中 碰撞 参数 所 决定 的 单 参数 轨道 族 进行 了 系统 的 数值 研究 ， 关 于 具体 的 计 
FAEN I ZZA 节 中 的 相关 介绍 。 这 族 单 参数 轨道 在 初始 时 刻 对 应 着 两 个 圆 轨 道 ， 而 碰撞 
参数 h 即 为 两 个 卫星 的 圆 轨 道 的 初始 轨道 半 长 径 之 差 。 当 碰撞 参数 hh 较 小 时 ， 数 值 计 算 发 
By 相对 运动 的 轨迹 是 一 个 几乎 对 称 的 马蹄 型 轨道 (Janus 和 Epimetheus 的 交换 轨道 的 
稳定 性 源 于 系统 的 h 值 很 小 )。 同 时 ， 假 设 h 为 小 量 ， 从 Hill 方程 可 简单 得 到 一 个 U 形 轨 迹 
的 代数 方程 ””， 由 该 方程 可 简单 解 得 两 个 卫星 相对 距离 最 小 值 与 初始 碰撞 参数 的 关系 。 
根据 Hill 方程 的 雅 可 比 积 分 ， 若 两 个 卫星 初始 轨道 均 为 圆 轨道 ， 则 当 它 们 相互 作用 又 远离 
后 ， 两 个 卫星 的 轨道 半 长 径 之 差 一 定 大 于 或 等 于 初始 碰撞 参数 ， 轨 道 半 长 径 之 差 增 大 对 应 
着 卫星 的 相对 运动 的 势能 降低 ， 考 虑 Hill 方程 的 雅 可 比 积 分 守恒 ， 这 部 分 减 小 的 能 量 转化 
PI 为 轨道 偏心 率 的 增 大 ， 即 相对 运动 的 动能 增 大 。 换 言 之 ， 由 碰撞 参数 单 参数 决定 的 两 个 初 
- 始 圆 轨道 的 相对 运动 可 分 为 两 种 情况 : (1) 对 称 的 马蹄 型 轨道 ， 且 相互 作用 后 轨道 偏心 率 

230; (2) 相互 作用 后 轨道 半 长 径 之 差 增 大 ， 且 轨道 偏心 率 不 为 零 。 事 实 上 ， 人 情况 (1) 对 应 
区 2 可 节 中 提 到 的 一 个 结论 4h BNI ( 记 < ho)， 相 对 运动 的 约 化 偏心 率 er 是 系统 的 
一 个 绝热 不 变量 。 当 初始 hh 较 大 时 (id h > hi)， 两 个 卫星 之 间 的 引力 始终 是 个 小 量 ， 此 时 
可 用 摄 动 方法 计算 两 个 轨道 在 交汇 前 后 的 微小 变化 。 当 九 在 上 述 ho A hs 中 间 变 化 时 
(BI ho <h<hi)， 轨 道 的 相对 运动 表现 出 十 分 复杂 的 演化 和 似乎 难以 穷尽 的 细节 。 具 体 而 
言 ， 当 在 某 个 临界 值 附 近 时 ， 相 对 于 h 的 微小 变化 ， 系 统 的 轨迹 却 表 现 出 不 连续 的 “ 突 
AE": 两 个 卫星 在 相互 作用 后 再 次 分 开 的 运动 方向 会 发 生 突变 ， 而 这 些 “ 突 变 ” 源 于 系统 中 
的 不 稳定 周期 轨道 。 具 体 地 说 ，h 的 某 个 临界 值 对 应 的 相对 运动 轨道 恰好 处 在 系统 某 个 周期 
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轨道 的 不 稳定 流 形 上 ， 当 疡 在 该 临界 值 附近 变 化 时 ， 轨 道 会 从 不 稳定 流 形 的 两 侧 渐 近 逼近 
该 周期 轨道 ， 进 而 从 不 同方 向 离开 ， 这 导致 轨迹 的 运动 方向 发 生 了 “突变 ” 系统 实际 的 运 
动 比 上 述 的 物理 图 像 更 为 复杂 ， 当 轨道 从 一 个 周期 轨道 的 不 稳定 流 形 离 开 后 ， 可 能 会 再 次 进 
入 另 一 个 周期 轨道 的 不 变 流 形 ; 对 于 轨道 族 的 演化 而 言 ， 已 有 的 计算 表明 在 一 个 有 限 的 
区 间 内 的 全 部 临界 值 似乎 是 不 可 数 的 零 测 集 ， 且 表现 出 复杂 的 自 相似 结构 ™。 因 为 h 的 临 
界 值 对 应 的 相对 运动 处 在 周期 轨道 的 渐 近 逼近 轨道 上 ， 这 对 应 着 系统 的 永久 俘获 ， 而 由 于 到 
的 临界 值 是 个 零 测 集 ， 故 永久 俘获 实际 上 几乎 不 会 发 生 。 

下 面 对 共 轨 运 动 的 长 期 演化 理论 做 简要 介绍 。 在 椭圆 型 限制 性 三 体 问题 下 ，érdi™™ 对 
三 角 平 动 点 附近 的 运动 构造 了 一 套 分 析 解 ， 基 于 该 解 6rdi™”™ 给 出 了 三 角 平 动 点 附近 运动 
KPRI EZ. Morais” 汪 用 经 典 方法 对 摄 动 函数 展开 到 二 阶 ( 见 节 )， 并 沿 着 
1/1 共振 的 可 积 近 似 的 运动 轨迹 对 哈密 顿 函数 平均 化 ， 进 而 得 到 1/1 共振 的 长 期 演化 方程 。 
斌 断 至 小 参数 的 二 阶 量 时 ， 平 均 化 运动 方程 可 以 显 式 地 解 出 ， 该 解 中 的 频率 仅 依赖 于 一 个 
表征 无 摄 1/1 共振 模型 中 轨迹 振幅 的 参数 ， 而 不 依赖 于 偏心 率 和 轨道 倾角 ， 故 Morais 称 之 
为 线性 长 期 理论 。Morais 进一步 研究 了 中 心 天 体 的 扁 率 以 及 第 三 体 引力 的 平均 化 效应 ， 并 
给 出 了 相应 的 摄 动 解 ， 该 解 包含 长 期 共振 对 应 的 小 分 母 。 应 用 该 线性 长 期 摄 动 理论 ，Morais 
计算 了 土星 摄 动 下 木星 Trojan 小 行星 以 及 天 王 星 的 卫星 系统 中 长 期 共振 的 位 置 ， 并 与 数值 
只 分 的 结果 进行 了 对 比 。 可 以 看 出 ，Morais 的 线性 长 期 摄 动 理论 适用 于 小 偏心 率 小 轨道 倾 
角 的 运动 ， 同 时 平均 化 的 有 效 性 要 求 了 1/1 共振 的 规则 运动 (远离 不 同 共振 态 之 间 的 边界 )， 
同时 1/1 共振 周期 要 远 小 于 长 期 演化 的 特征 时 标 。 基 于 小 偏心 率 e 及 轨道 倾角 i 的 经 典 摄 动 
函数 展开 ”，Robutel 和 Pousse” 对 行星 三 体 问题 的 共 轨 运动 提出 了 一 个 理想 共振 模型 ， 与 
Yoder 和 Morais 等 人 模型 的 主要 区 别 在 于 该 模型 具有 五 个 平衡 点 。Robutel 和 Pousse 进 
步 研究 理想 共振 模型 中 不 同 轨迹 的 长 期 演化 方程 类似 于 Morais 的 线性 长 期 理论 ， 截 断 至 
O(e?) 时 该 方程 是 系数 为 周期 函数 的 线性 微分 方程 组 。 由 于 考虑 小 偏心 率 和 小 轨道 倾角 ， 长 
期 演化 方程 此 时 等 价 于 理想 共振 模型 不 同方 向 上 的 变 分 方程 。 在 平衡 点 处 ， 该 方程 组 可 直 
接 解 出 ， 在 理想 共振 模型 的 任 一 周期 轨道 上 ， 方 程 组 可 根据 Floquet 理论 解 出 。 方 程 组 解 的 
频率 对 应 着 1/1 共振 轨迹 的 长 期 演化 频率 ， 由 此 Robutel 和 Pousse 计算 了 1/1 共振 中 可 能 
发 生 的 次 级 共振 的 位 置 。 特 别 地 ， 在 点 LJ4“ 偏 心率 方向 ”的 变 分 方程 有 一 个 零 特 征 根 和 一 
= 个 非 零 特征 根 ， 零 特征 根 对 应 着 拉 格 朗 日 周期 轨道 ， 非 零 特 征 根 对 应 着 所 谓 的 反 拉 格 朗 日 
= 周期 轨道 ”"， 在 点 L 有 类 似 的 结果 。 需 要 指出 的 是 ， 这 些 周 期 轨道 是 定义 在 平均 化 系统 中 。 
FH Lyapunov 中 心 定 理 ， 点 Za 处 生发 了 一 族 与 零 特 征 根 的 特征 方向 相 切 的 周期 轨道 ， 即 拉 
格 朗 日 周期 轨道 。 但 由 于 另 一 特征 根 为 零 ，Lyapunov 中 心 定 理 无 法 保证 反 拉 格 朗 日 轨道 族 
的 存在 性 。Robutel 和 Pousse 进而 对 点 Ly 附近 的 运动 构造 了 高 阶 的 Birkhoff normal form, 
并 证 明了 从 点 La 生发 了 一 族 反 拉 格 朗 日 周期 轨道 。 

上 面 介绍 的 长 期 演化 理论 仅 适用 于 近 贺 近 平面 的 共 轨 运动 。Namouni” 在 研究 中 首先 指 
出 : 在 大 偏心 率 或 高 轨道 倾角 处 存在 着 与 熟知 的 1/1 共振 轨迹 不 同 的 共振 态 ， 例 如 蔬 Z 避 市 
中 介绍 的 跃迁 (transition) 轨道 等 。Namouni 首先 在 Hill 三 体 问 题 中 发 现 了 轨道 因 长 期 演化 
在 不 同 的 1/1 共振 态 之 间 迁 移 的 可 能 性 ， 进 而 用 数值 积分 方法 研究 了 一 些 选 定 轨 道 的 长 期 
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演化 。 在 相 图 上 的 不 同 区 域 ，Namouni 确认 了 与 轨道 长 期 演化 相关 的 动力 学 现象 ， 例 如 围 
SE wr = 0° 或 180。 的 Kozai 共振 ， 并 发 现 这 些 束缚 在 长 期 共振 的 轨道 往往 对 应 某 一 类 跃迁 
轨道 。Namouni 等 人 指出: 可 能 存在 着 尚未 发 现 的 运行 在 1/1 共振 跃迁 轨道 上 的 近 地 小 
行星 ， 这 些小 行星 有 可 能 被 地 球 短暂 地 俘获 在 拟 卫 星 轨 道上 ”。 同 时 ，Namouni 等 人 指出 ， 
不 同 1/1 共振 态 之 间 的 跃迁 也 为 地 球 防御 近 地 小 行星 的 撞击 提供 了 一 种 保护 机 制 。 通 过 将 摄 
动 函数 在 非 零 偏心 率 和 轨道 倾角 处 展开 以 及 数值 平均 ，Nesvorny 等 人 "用 哈密 顿 系统 研究 
了 两 种 情形 的 共 轨 运动 : 平面 椭圆 型 限制 性 三 体 问题 及 空间 圆 型 限制 性 三 体 问题 ， 两 个 问 
题 均 为 三 自由 度 。 通 过 数值 平均 (UELLE) 去 除 哈 密 顿 函数 中 的 快 变量 ， 系 统 降 为 二 自由 
度 ; 对 平均 化 系统 应 用 绝热 摄 动 理论 ( 见 区 [L 引 节 )， 系 统 降 为 一 自由 度 ， 对 应 原 系统 的 长 期 
演化 。 通 过 将 二 自由 度 平 均 化 系统 中 的 慢 变 量 “ 冻 结 ”，Nesvorny 等 人 给 出 了 相 空 间 不 同 区 
域 的 1/1 共振 相 图 ， 并 发 现 了 不 同 于 平面 情形 的 1/1 共振 态 (例如 图 苞 ))， 并 指出 了 在 轨道 
的 长 期 进 动 下 ， 系 统 在 不 同 1/1 共振 态 之 间 迁 移 的 可 能 性 ; 一 自由 度 系 统 给 出 的 长 期 演化 相 
图 给 出 了 与 Namouni” 工作 中 相近 的 动力 学 性 质 ，Nesvorny 等 人 用 数值 积分 对 长 期 演化 相 
图 上 的 结构 进行 了 验证 ， 并 指出 在 不 同 共振 态 间 迁移 的 轨道 基本 是 混沌 的 。 
Barrabés 和 Mikkola” 计算 了 限制 性 问题 中 平面 和 空间 中 的 对 称 马 蹄 型 周期 轨道 : 周 
期 轨道 的 迭代 与 延 拓 使 用 了 Gómez fll Mondelo 给 出 的 一 般 方法 ， 而 周期 轨道 的 初始 值 出 
Llibre 和 Ollé” 的 方法 给 出 。 通 过 计算 ，Barrabés 和 Mikkola 研究 了 平面 马蹄 型 周期 轨道 
族 的 结构 和 演化 ， 发 现 了 轨道 倾角 i = 17° 的 稳定 的 周期 轨道 。Barrabés 和 Ol16” 系统 地 研 
究 了 平面 限制 性 三 体 问题 中 的 马蹄 型 周期 轨道 族 。 马 蹄 型 轨道 在 二 体 问 题 极 限 下 (ú 一 0) 对 
应 着 两 条 开 普 勒 轨道 ， 对 称 马 蹄 型 周期 轨道 要 求 这 两 条 开 普 勒 轨道 在 转动 坐标 系 下 也 为 对 
称 周期 轨道 。Barrabés 和 Ollé 导出 了 这 族 开 普 勒 轨道 的 初始 条 件 ， 并 称 之 为 马蹄 型 轨道 的 
生成 族 。 基 于 上 述 的 生成 族 ，Barrabés 和 Olé 通过 数值 延 拓 得 到 了 u > 0 Ap 32 s (不 妨 
设 该 值 为 = uu) 的 马蹄 型 周期 轨道 族 ， 并 根据 生成 族 的 结构 解释 了 = He 时 马蹄 型 周期 
轨道 族 的 特征 。 在 得 到 j= us 处 Ls 点 附近 的 周期 轨道 族 后 ， 一 个 自然 的 想法 便 是 以 这 族 
轨道 为 初始 值 ， 进 一 步 将 周期 轨道 族 延 拓 至 u = pa 十 6， 其 中 6 充分 小 使 得 迭代 收敛 。 然 
而 ，Barrabés 和 Ollé 通过 计算 发 现 L, 点 附近 的 不 变 流 形 非常 敏感 地 依赖 于 TEL IU] NAE 
化 ， 这 使 得 简单 地 沿 着 质量 参数 u 进行 延 拓 变 得 困难 。 根 据 上 述 特点 ， 为 分 析 Vu € (0,1/2] 
= 系统 中 点 Ls 附近 的 HPO, Barrabés 和 Ollé 根据 以 下 机 制 系统 地 研究 了 点 La 的 不 稳定 流 
w 的 形状 当 点 Ls 的 不 稳定 流 形 为 马蹄 型 的 同 宿 轨道 时 ( 记 此 时 /= jwwomo)， 在 该 同 宿 轨 
道 的 邻 域内 可 找到 无 穷 多 个 马蹄 型 周期 轨道 。 通 过 计算 ，Barrabés 和 Olé 发 现 了 两 个 极限 
为 零 的 aome 的 无 穷 序 列 ， 而 对 Uomo 值 附近 的 系统 的 点 Ls 不 变 流 形 的 计算 显示 : 系统 存 
在 双 “ 周 期 ”的 杞 蹄 型 同 宿 轨道 ， 即 该 同 宿 轨道 在 一 个 “周期 ”内 绕 转 点 L. 和 Ls 两 次 。 
Barrabés 和 Ollé 根据 以 上 机 制 解释 了 Vu e (0,1/2] 系统 中 点 Ls 附近 的 HPO 的 存在 性 。 当 
然 ， 系 统 存在 其 他 机 制 “ 生 成 ”马蹄 型 周期 轨道 ， 如 Deprit 和 Henrard"" 的 计算 显示 ， 点 
Ls 的 不 稳定 周期 轨道 的 不 变 流 形 存在 马蹄 型 同 宿 轨 道 。Barrabks 和 Ollé 发 现 : 点 La 的 不 
稳定 流 形 在 某 些 j, 值 处 会 十 分 靠近 点 L 和 LK。， 甚 至 存在 与 第 二 主 天 体 接触 的 碰撞 流 形 ， 
这 表明 马蹄 型 轨道 的 演化 与 系统 全 局 的 动力 学 结构 密切 相关 。Hou 和 Liu” 通过 数值 研究 
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HH: 对 称 HPO 可 由 围绕 点 Ls 的 平面 Lyapunov 周期 轨道 族 的 m 分 又 轨道 得 到 ， 即 对 称 
HPO 是 连接 这 些 分 叉 周期 轨道 的 周期 轨道 族 。 

Hadjidemetriou 等 人 指出 平面 行星 三 体 问题 共 轨 运动 中 存在 两 族 对 称 周期 轨道 (一 
族 稳定 ， 一 族 不 稳定 )， 这 两 族 轨道 在 分 析 上 可 由 点 Ls 处 的 变 分 方程 得 到 "。 在 稳定 的 周 
期 轨道 上 ， 两 个 行星 的 轨道 偏心 率 初始 时 近似 满足 me = m2e2””， 近 点 角 距 相差 180° 
与 中 心 天 体 共 线 ， 初 始 平 经 度 也 相差 180?， 这 一 构 型 对 应 着 拟 卫星 轨道 ”。 在 不 稳定 的 
周期 轨道 上 ， 两 个 行星 的 轨道 偏心 率 满足 e = ez， 近 点 角 距 相差 0° ( 即 在 中 心 天 体 同一 
fl) 且 与 中 心 天 体 共 线 ， 初 始 平 经 度 也 相差 180"， 这 一 构 型 对 应 着 Euler 的 共 线 平 动 点 。 
Hadjidemetriou 等 人 " 用 庞 加 莱 截 面 计算 了 周期 轨道 族 邻 域内 的 运动 。 从 2/1 和 3/1 平 运动 
t 振 的 研究 中 己 经 知道 ”在 非 保守 力作 用 下 行星 系统 可 以 沿 着 稳定 的 周期 轨道 族 缓慢 
演化 。Hadjidemetriou 和 Voyatzis” 发 现 ， 当 一 个 拖 忠 力作 用 在 1/1 共振 系统 中 时 ， 系 统 可 
沿 着 上 述 的 稳定 周期 轨道 族 缓慢 地 由 行星 轨道 向 卫星 轨道 迁移 ， 这 提供 了 一 个 由 行星 系统 
产生 卫星 系统 的 机 制 。 

2017 年 ，Wiegert 5& A^ ”证 实 太 阳 系 中 存在 一 颗 与 木星 处 于 1/1 逆行 共振 的 小 行星 
2015 BZ509。 从 会 合 坐标 系 下 看 ， 小 行星 2015 BZ509 的 轨道 形状 为 出 形 ， 数 值 计 算 表 明 
2015 BZ509 的 轨道 在 百 万 年 尺度 上 是 稳定 的 。 对 于 逆行 1/1 共振 的 研究 开始 于 2015 BZ509 
的 发 现 之 前 。Morais 和 Namouni" 对 1/1 逆行 共振 建立 了 分 析 和 半分 析 模型 ， 并 发 现 了 几 
种 不 同 的 共振 态 。Morais 和 Namouni” 通 过 数值 计算 轨道 稳定 性 指标 发 现 ， 在 圆 型 限制 性 
问题 下 存在 稳定 的 空间 1/1 道行 共振 轨道 。Huang 等 人 "研究 了 1/1 逆行 共振 中 的 Kozai 
t 振 。 在 空间 圆 型 限制 性 问题 下 ，Sidorenko ”用 哈密 顿 系统 构造 了 1/1 逆行 共振 的 长 期 摄 
动 理 论 。 


= 
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本 文 回顾 了 限制 性 及 行星 三 体 问 题 中 的 共 轨 运 动 相 关 的 模型 、 方 法 及 动力 学 研究 。 第 中 
章 在 物理 直观 上 解释 了 共 轨 运动 及 1/1 平 运动 共振 的 含义 ， 介 绍 了 限制 性 及 行星 三 体 问题 
中 典型 的 共 轨 运动 类 型 ， 并 对 太阳 系 中 的 自然 共 轨 系统 、 太 阳 系 外 共 轨 系统 的 探测 现状 做 
了 简要 介绍 。 第 加 章 主要 介绍 了 处 理 共 轨 运动 的 两 种 常用 方法 : 哈密 顿 系统 以 及 Hull 三 体 
问题 。 哈密 顿 系 统 可 以 研究 1/1 共振 全 局 或 局 部 的 动力 学 性 质 ， 区 II 节 给 出 了 限制 性 及 行 
星 三 体 问 题 的 哈密 顿 函 数 ， 区 TI3 节 介绍 了 适用 于 1/1 共振 的 摄 动 函数 展开 方法 ， 并 给 出 了 
经 典 展开 方法 截断 至 二 阶 量 的 结果 ; 区 LI3 节 介绍 了 哈密 顿 系统 降低 自由 度 的 两 种 方法 : 1E 
规 化 和 平均 化 ， 并 给 出 相关 算 例 ,区 工 了 节 以 行星 三 体 问题 为 例 介 绍 了 1⁄1 共振 的 理想 共振 
模型 (又 称 为 可 积 近似 ); 区 TI 引 节 首 先 介绍 了 经 典 的 3/1 共振 长 期 摄 动 理论 (主要 是 Wisdom 
及 Henrard 等 人 的 工作 )， 进 而 讨论 了 1/1 共振 的 长 期 摄 动 理 论 。Hill 的 三 体 问 题 是 处 理 共 
轨 和 运动 的 另 一 有 效 途 径 ， 区 2 节 给 出 了 经 典 的 Hill 方程 ， 并 推导 了 适用 行星 三 体 问题 的 更 
一 般 的 椭圆 型 Hill 方程 ， 苞 Z 习 节 给 出 了 圆 型 Hill 方程 的 无 摄 解 、U 型 轨迹 分 析 解 ， 并 介绍 
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了 连接 Hill 方程 外 部 近似 与 内 部 近似 的 四 参数 的 渐进 逼近 解 ; 区 Z 引 节 推 导 了 圆 型 Hill 方程 
的 正则 形式 ， 为 Hill 方 程 的 平均 化 做 了 准备 ， 基 于 圆 型 Hill 方程 的 正则 形式 ， 区 2 可 节 介 绍 了 
Namouni 处 理 共 轨 运 动 的 摄 动 方法 。 最 后 ， 区 3 引 节 简单 介绍 了 基于 轨道 数值 积分 的 常用 的 数 
值 方法 : 周期 轨道 延 拓 、 庞 加 莱 截 面 、 稳 定性 指标 的 计算 等 。 第 回 章 回顾 了 共 轨 运动 相关 的 
动力 学 研究 ， 包 括 晴 蛙 型 轨道 演化 的 Brown 猜想 ， 行 星 三 体 问题 的 1/1 理想 共振 模型 ,“ 卫 
星 遭 遇 ” 问 题 中 的 数值 计算 ， 小 偏心 率 及 轨道 倾角 下 的 共 轨 运动 的 线性 长 期 理论 ， 大 偏心 率 
和 大 轨道 倾角 下 新 的 1/1 共振 态 以 及 不 同 共振 态 之 间 的 跃迁 ， 限 制 性 及 行星 三 体 问题 中 周期 
轨道 族 和 不 变 流 形 的 计算 等 。 
第 回 章 提 到 ， 哈 密 顿 系统 仍 有 一 些 技术 性 问题 尚未 解决 ， 其 中 之 一 是 摄 动 函数 展开 在 靠 
近 磁 撞 奇 点 时 的 收敛 性 问题 。Pousse 等 人 ” 指出， 平均 化 方法 在 拟 卫 星 区 域 的 有 效 性 受到 限 
制 的 原因 有 两 个 ， 一 是 限制 性 三 体 问题 中 与 拟 卫 星 轨道 对 应 的 / 族 周期 轨道 "在 第 二 主 天 体 
附近 的 频率 并 不 存在 短 周 期 与 长 周期 之 分 ， 故 平均 化 方法 从 根本 上 不 适用 ; 第 二 个 原因 是 
技术 性 的 ， 即 摄 动 函数 展开 在 靠近 碰撞 奇 点 处 的 收敛 问题 ， 而 这 一 收敛 性 问题 可 能 通过 新 
的 展开 方法 得 到 改善 。 哈 密 顿 系统 的 另 一 技术 性 问题 是 摄 动 函 数 展 开 受 限于 小 偏心 率 及 小 
轨道 倾角 ，Brumberg 等 人 "提出 了 用 椭圆 函数 描述 开 普 勒 轨道 的 方法 ， 这 为 非 零 偏 心率 
处 摄 动 函 数 展开 提供 了 一 个 分 析 途 径 。 区 工 可 节 提 到 ，Garfinkel 通过 他 的 广义 理想 共振 模型 
对 限制 性 问题 的 共 轨 运动 得 到 了 二 阶 解 ， 用 相同 的 方法 可 以 尝试 处 理 行星 三 体 问题 。Henon 
和 Petit^" 以 及 Spirig 和 Waldvogel” 给 出 了 平面 圆 型 Hill 方程 的 一 族 四 参数 渐 近 逼近 解 ， 
可 以 考虑 将 该 解 推广 至 空间 情形 下 或 椭圆 型 Hill 方程 。 从 数值 计算 方法 出 发 ， 可 进一步 考 
虑 平面 圆 型 Hill 方程 下 初始 轨道 为 椭圆 轨道 的 “卫星 遭遇 ”问题 ”， 以 及 初始 轨道 空间 Z 
方向 分 量 不 为 零 的 情形 。 在 平面 加 型 限制 性 问题 中 ， 由 平均 化 分 析 理 论 得 到 的 运动 性 态 CF 
衡 点 、 周 期 轨道 ) 与 完整 系统 中 数值 方法 得 到 的 周期 轨道 族 、 拟 周期 轨道 之 间 的 联系 是 清楚 
的 所。 类 似 地 ， 可 以 考虑 大 偏心 率 或 高 倾角 处 以 及 非 限 制 性 问题 中 ， 周 期 轨道 族 以 及 拟 周期 
轨道 与 分 析 理 论 之 间 的 联系 。 
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表达 式 分 别 为 : 
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其 中 27 一 1 
p = (ai + as = 20102 COS (Ai = à2)) 2 
参考 文献 
1] Hénon M, Guyot M. In: Giacaglia G E O, ed. Periodic Orbits Stability and Resonances. Dordrecht: 
Springer, 1970: 349 
2] Pousse A, Robutel P, Vienne A. Celestial Mechanics and Dynamical Astronomy, 2017, 128(4): 383 
3] Nesvorny D, Thomas F, Ferraz-Mello S, et al. Celestial Mechanics and Dynamical Astronomy, 2002, 
82(4): 323 
4] Namouni F. Icarus, 1999, 137(2): 293 
5] Laskar J, Robutel P. Celestial Mechanics and Dynamical Astronomy, 1995, 62(3): 193 
6] Robutel P, Pousse A. Celestial Mechanics and Dynamical Astronomy, 2013, 117(1): 17 
7| Hadjidemetriou J D, Psychoyos D, Voyatzis G. Celestial Mechanics and Dynamical Astronomy, 2009, 
104(1-2): 23 
8] Giuppone C A, Beaugé C, Michtchenko T A, et al. MNRAS, 2010, 407(1): 390 
9] Greenstreet S, Gladman B, Ngo H. AJ, 2020, 160(3): 144 
10] Laughlin G, Chambers J E. AJ, 2002, 124(1): 592 
11] Giuppone C A, Benitez-Llambay P, Beaugé C. MNRAS, 2012, 421(1): 356 
12] GoZdziewski K, Konacki M. ApJ, 2006, 647(1): 573 
13] Rein H. MNRAS, 2015, 448: L58 
14] Lissauer J J, Dawson RI, Tremaine S. Nature, 2014, 513(7518): 336 
15] Perryman M. The Exoplanet Handbook. Cambridge: Cambridge University Press, 2018: 76 
16] Leleu A, Coleman G A L, Ataiee S. A&A, 2019, 631: A6 
17] Ferraz-Mello S. Canonical Perturbation Theories - Degenerate Systems and Resonance. Vol. 345. New 
York: Springer, 2007: 174 
18] Celletti A, Chierchia L. Celestial Mechanics and Dynamical Astronomy, 2006, 95(1-4): 117 
19] Kaula WM. AJ, 1962, 67: 300 
20] Harrington RS. AJ, 1968, 73: 190 
21] Brouwer D, Clemence G M. Methods of Celestial Mechanics. New York: Academic Press Inc, 1961: 60 
22] Morais MH M. A&A, 1999, 350: 318 
23] oW. 卫星 轨道 力学 算法 . 南京 : 南京 大 学 出 版 社 , 2019 
24] Kaula W M. Geophysical Journal, 1961, 5(2): 104 
25] Robutel P, Niederman L, Pousse A. arXiv e-prints, 2015: arXiv:1506.02870 
26] Giuppone C A, Leiva AM. MNRAS, 2016, 460(1): 966 
27] Ferraz-Mello S. A&A, 1987, 183(2): 397 
28] Ferraz-Mello S, Sato M. A&A, 1989, 225: 541 
29] Roig F, Simula A, Ferraz-Mello S, et al. A&A, 1998, 329: 339 
30] Páez RI, Locatelli U, Efthymiopoulos C. Celestial Mechanics and Dynamical Astronomy, 2016, 126(4): 
519 
31] Jorba A, Masdemont J. Physica D Nonlinear Phenomena, 1999, 132(1-2): 189 
32] Tan P, Hou X Y, Liao X H. Celestial Mechanics and Dynamical Astronomy, 2021, 133(7): 31 
33] Arnold V I. Mathematical methods of classical mechanics. New York: Springer, 1978: 292 
34] Deprit A. Celestial Mechanics, 1969, 1(1): 12 
35] Garfinkel B. AJ, 1966, 71: 657 
36] Cors J M, Palacian J F, Yanguas P. SIAM Journal on Applied Dynamical Systems, 2019, 18(1): 334 


202306.00372v1 


" 
m 


chinaXiv 


1 期 


37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 


50 
51 
52 
53 
54 
55 
56 


57 
58 
59 
60 
61 


62 
63 


64 
65 
66 
67 
68 
69 
70 
val 
72 
73 


74 
75 
76 
77 


78 


ChinaXiv 合 作 期 


by, 5E. 限制 性 及 行星 三 体 问题 中 的 共 轨 运动 研究 回顾 93 


Niederman L, Pousse A, Robutel P. Communications in Mathematical Physics, 2020, 377(1): 551 
Garfinkel B. Celestial Mechanics, 1976, 13(2): 229 

Garfinkel B. AJ, 1977, 82: 368 

Garfinkel B. Celestial Mechanics, 1978, 18(3): 259 

Wisdom J. Icarus, 1985, 63(2): 272 

Henrard J. Celestial Mechanics and Dynamical Astronomy, 1990, 49(1): 43 

Murray C D, Dermott S F. Solar System Dynamics. Cambridge: Cambridge University Press, 1999: 462 
Lichtenberg A J, Lieberman M A. Regular and stochastic motion. New York: Springer, 1983: 88 
Henrard J, Lemaitre A. Celestial Mechanics, 1986, 39(3): 213 

Henrard J, Caranicolas N D. Celestial Mechanics and Dynamical Astronomy, 1990, 47(2): 99 

Beaugé C, Roig F. Icarus, 2001, 153(2): 391 

Kozai Y. Celestial Mechanics, 1985, 36(1): 47 

Giacaglia G E O, Nacozy P E. In: Giacaglia G E O, ed. Periodic Orbits Stability and Resonances. 
Dordrecht: Springer, 1970: 96 

Li M, Huang Y K, Gong S P. Research in Astronomy and Astrophysics (RAA), 2021, 21(2): 25 

Morais MH M. A&A, 2001, 369: 677 

Beauge C, Ferraz-Mello S. Icarus, 1994, 110(2): 239 

Lhotka C, Celletti A. Icarus, 2015, 250: 249 

Sicardy B, Dubois V. Celestial Mechanics and Dynamical Astronomy, 2003, 86(4): 321 

Henon M, Petit J M. Celestial Mechanics, 1986, 38(1): 67 

Spirig F, Waldvogel J. In: Szebehely V G, ed. Stability of the Solar System and its Minor Natural and 
Artificial Bodies. Dordrecht, Holland: D. Reidel Publishing Company, 1985: 53 

Petit J M, Henon M. Icarus, 1986, 66(3): 536 

Dermott S F, Murray C D. Icarus, 1981, 48(1): 1 

Kruskal M. Journal of Mathematical Physics, 1962, 3(4): 806 

Wisdom J. AJ, 1982, 87: 577 

Sandor Z, Morais MH M. In: Pretka-Ziomek H, Wnuk E, Seidelmann P K, et al, eds. Dynamics of 
Natural and Artificial Celestial Bodies. Dordrecht: Springer, 2002, 81: 263 

Hadjidemetriou J D. Predictability, Stability, and Chaos in N-Body Dynamical Systems. New York: 
Plenum Press, 1991: 157 

Henon M. Generating Families in the Restricted Three-Body Problem. Berlin, Heidelberg: Springer, 
1997: 1 

Rabe E. AJ, 1961, 66: 500 

Taylor D B. A&A, 1981, 103(2): 288 

Barrabés E, Ollé M. Nonlinearity, 2006, 19(9): 2065 

Kreisman B B. Cosmic Research, 2016, 54(6): 475 

Morais M H M, Namouni F. MNRAS, 2019, 490(3): 3799 

Morais M H M, Namouni F. Celestial Mechanics and Dynamical Astronomy, 2013, 117(4): 405 

Zhao Z Y, Liu L. Icarus, 1992, 100(1): 136 

Funk B, Dvorak R, Schwarz R. Celestial Mechanics and Dynamical Astronomy, 2013, 117: 1 

Morais M H M, Namouni F. Celestial Mechanics and Dynamical Astronomy, 2016, 125(1): 91 

Skokos C. The Lyapunov Characteristic Exponents and Their Computation. Berlin, Heidelberg: Springer, 
2010: 63 

Cincotta P M, Simó C. A&AS, 2000, 147: 205 

Laskar J. arXiv Mathematics e-prints, 2003: math/0305364 

Brown E W. MNRAS, 1911, 71: 438 

Deprit A, Henrard J. In: Giacaglia G E O, ed. Periodic Orbits Stability and Resonances. Dordrecht: 
Springer, 1970: 1 

Henrard J. Celestial Mechanics, 1983, 31(2): 115 


Ti 


94 天 文学 进展 40 卷 


79] Garfinkel B. Celestial Mechanics, 1985, 36(1): 19 

80] Aksnes K. In: Szebehely V G, ed. Stability of the Solar System and its Minor Natural and Artificial 
Bodies. Dordrecht, Holland: D. Reidel Publishing Company, 1985: 3 

81] Dermott S F, Murray C D. Icarus, 1981, 48(1): 12 

82] Yoder C F, Colombo G, Synnott S P, et al. Icarus, 1983, 53(3): 431 

83] Goldreich P, Tremaine S. ApJ, 1980, 241: 425 

84] Petit JM, Hénon M. On a Cantor structure in a satellite scattering problem. Vol. 355. Berlin, Heidelberg: 
Springer, 1990: 225 

85] Erdi B. Celestial Mechanics, 1977, 15(3): 367 

86] Erdi B. Celestial Mechanics, 1978, 18(2): 141 

87] Erdi B. Celestial Mechanics, 1979, 20(1): 59 

88] Erdi B. Celestial Mechanics, 1981, 24(4): 377 

89] Namouni F, Christou A A, Murray C D. Phys. Rev. Lett., 1999, 83(13): 2506 

90] Mikkola S, Innanen K, Wiegert P, et al. MNRAS, 2006, 369(1): 15 

91] Barrabés E, Mikkola S. A&A, 2005, 432(3): 1115 

92] Gómez G, Mondelo J M. Physica D Nonlinear Phenomena, 2001, 157(4): 283 

93] Llibre J, Ollé M. A&A, 2001, 378: 1087 

94] Deprit A, Henrard J. AJ, 1969, 74: 308 

95] Hou X Y, Liu L. AJ, 2008, 136(1): 67 

96] Ferraz-Mello S, Beaugé C, Michtchenko T A. Celestial Mechanics and Dynamical Astronomy, 2003, 87(1): 
99 

97] Beaugé C, Michtchenko T A, Ferraz-Mello S. MNRAS, 2006, 365(4): 1160 

98] Hadjidemetriou J D, Voyatzis G. Celestial Mechanics and Dynamical Astronomy, 2011, 111(1-2): 179 
99] Wiegert P, Connors M, Veillet C. Nature, 2017, 543(7647): 687 

[100] Morais H, Namouni F. Nature, 2017, 543(7647): 635 

[101] Huang Y, Li M, Li J, et al. MNRAS, 2018, 481(4): 5401 

[102] Sidorenko V V. AJ, 2020, 160(6): 257 

[103] Brumberg E, Brumberg V A, Konrad T, et al. Celestial Mechanics and Dynamical Astronomy, 1995, 
61(4): 369 


A Review on Co-orbital Motion in Restricted and Planetary 
Three-body Problem 
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Abstract: The 1/1 mean motion resonance may be referred to as the lowest order mean 
motion resonance in the restricted or planetary three-body problem. The well-known five 
libration points in the circular restricted three-body problem correspond to the five equilib- 
riums of the 1/1 resonance. Trapped in different 1/1 libration modes, the system moves on 


different types of trajectory in space, such as the tadpole orbits and the horseshoe orbits 
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in case of small orbital eccentricity and inclination. In case of moderate or large eccentric- 
ity and inclination, it is known that other novel 1/1 libration modes also exist. Coorbital 
bodies are not rare in the solar system, for example the Trojans asteroids and the coorbital 
satellite systems of Saturn. Recently, dozens of coorbital bodies with Earth as the host have 
been identified among the near-Earth asteroids. Due to the orbital precession, the planetary 
graviational perturbations and other possible effects, these coorbital bodies are believed to 
transit recurrently between different 1/1 libration modes before their eventual escapes. The 
Hamiltonian system and the Hill's three-body problem are two effective approaches in the 
studies of coorbital motion. To apply the perturbation theory to the Hamiltonian system, 
one first develops the disturbing function into an appropriate form, then uses the normal- 
ization procedure or the averaging technique to eliminate the short-period terms in order to 
lower the degree of freedom of the system. Then, one studies the integrable approximation 
(or referred to as the ideal resonance model), based on which the theory on the secular evo- 
lution may be constructed. The approach of Hamiltonian system can investigate the global 
dynamics of coorbital motion with suitable expansion of the disturbing function. On the 
other hand, the Hill's three-body problem focuses mainly on relative motion of two bodies 
when they are close. The Hill's equation in the context of the circular restricted three-body 
problem is well known. However, the general Hill's three-body problem whose equation of 
motion has exactly the same form can handle the non-restricted case where the mass of each 
body is non-negligible. Under certain assumptions, the equation of motion of the general 
Hill's three-body problem can be written in a “canonical” form. Based on this canonical 
form, the averaging technique can be used to construct a perturbation theory on the secular 
evolution. Besides the two analytical approaches, the numerical methods for the study of 
coorbital motion involve the continuation of periodic orbit, the surface of section and the 


calculation of invariant manifold or invariant tori around the equilibrium. 


Key words: three-body problem; 1/1 mean motion resonance; coorbital motion 


